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The Beginning and Evolution of the Universe 

Bharat Ratra 1 and Michael S. Vogeley 2 

ABSTRACT 

We review the current standard model for the evolution of the Universe from an early 
inflationary epoch to the complex hierarchy of structure seen today. We summarize and 
provide key references for the following topics: observations of the expanding Universe; 
the hot early Universe and nucleosynthesis; theory and observations of the cosmic mi- 
crowave background; Big Bang cosmology; inflation; dark matter and dark energy; 
theory of structure formation; the cold dark matter model; galaxy formation; cosmolog- 
ical simulations; observations of galaxies, clusters, and quasars; statistical measures of 
large-scale structure; and measurement of cosmological parameters. We conclude with 
discussion of some open questions in cosmology. This review is designed to provide a 
graduate student or other new worker in the field an introduction to the cosmological 
literature. 

Subject headings: Review 



1. Introduction 

It is the current opinion of many physicists that the Universe is well described by what Fred 
Hoyle termed a Big Bang Model, in which the Universe expanded from a denser hotter childhood to 
its current adolescence, with a present energy budget dominated by dark energy and less so by dark 
matter, neither of which have been detected in the laboratory, with the stuff biological systems, 
planets, stars, and all visible matter are made of (called baryonic matter by cosmologists) being a 
very small tracer on this dark sea, and with electromagnetic radiation being an even less significant 
contributor. Galaxies and groups and clusters of galaxies are locally distributed inhomogeneously 
in space, but on large enough scales and in a statistical sense the distribution approaches isotropy. 
This is supported by other electromagnetic distributions such as the X-ray and cosmic microwave 
backgrounds, which are close to isotropic. As one looks out further into space, as a consequence of 
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the finite speed of light, one sees objects as they were at earlier times, and there is clear observational 
evidence for temporal evolution in the distribution of various objects such as galaxies. 

At earlier times the Universe was hotter and denser, at some stage so hot that atoms could 
not exist. Nuclear physics reactions between protons, neutrons, etc., in the cooling expanding 
Universe resulted in the (nucleo) synthesis of the lighter elements (nuclei) such as D, 4 He, and 7 Li, 
with abundances in good accord with what is observed, and with the photons left over forming a 
residual cosmic microwave background (CMB) also in good agreement with what is observed. 

Given initial inhomogeneities in the mass distribution at an earlier time, processing of these 
by the expansion of the Universe, gravitational instability, pressure gradients, and microphysical 
processes, gives rise to observed anisotropics in the CMB and the current large-scale distribution 
of nonrelativistic matter; the situation on smaller spatial scales, where galaxies form, is murkier. 
Observations indicate that the needed initial inhomogeneities are most likely of the special form 
known as scale invariant, and that the simplest best-fitting Big Bang Model has flat spatial ge- 
ometry. These facts could be the consequence of a simple inflationary infancy of the Universe, a 
very early period of extremely rapid expansion, which stretched zero-point quantum-mechanical 
fluctuations to larger length scales and transmuted them into the needed classical inhomogeneities 
in the mass-energy distribution. At the end of the inflationary expansion all radiation and matter 
is generated as the Universe moves into the usual Big Bang Model epoch. Inflation has roots in 
models of very high-energy physics. Because of electromagnetic charge screening, gravity is the 
dominant large-scale force. General relativity is the best theory of gravity. 

This review attempts to elaborate on this picture. Given the Tantalus principle of cosmol- 
ogy (and most of astrophysics), that one can see but not "touch" — which makes this a unique 
field of physics — there have been many false starts and even much confusion and many missed 
opportunities along what most now feel is the right track. Given space constraints we cannot do 
justice to what are now felt to be false starts, nor will we discuss more than one or two examples 
of confusion and missed opportunities. We attempt here to simply describe what is now thought 
to be a reasonable standard model of cosmology and trace the development of what are now felt to 
be the important threads in this tapestry; time will tell whether our use of "reasonable standard" 
is more than just youthful arrogance (or possibly middle-aged complacence). 

In the following sections we review the current standard model of cosmology, with emphasis 
in parts on some historical roots, citing historically significant and more modern papers as well 
as review articles. We begin with discussion of the foundations of the Big Bang Model in Sec. 
[21 which summarizes research in the half century from Einstein's foundational paper on modern 
cosmology until the late 1960's discovery of the CMB radiation, as well as some loose ends. Section 
[3] discusses inflation, which provides an explanation of the Big Bang that is widely felt to be 
reasonable. Dark energy and dark matter, the two (as yet not directly detected) main components 
of the energy budget of the present Universe are reviewed in Sec. HI Further topics include the 
growth of structure in the Universe (Sec. [5]), observations of large-scale structure in the Universe 
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(Sec.[6j), and estimates of cosmological parameters (Sec. [7]). We conclude in Sec.[8]with a discussion 
of what are now thought to be relevant open questions and directions in which the field appears to 
be moving. 

We use hardly any mathematical equations in this review. In some cases this results in dis- 
guising the true technical complexity of the issues we discuss. 

We exclude from this review a number of theoretical topics: quantum cosmology, the multiverse 
scenario, string gas cosmology, braneworld and higher dimensional scenarios, and other modifica- 
tions of the Einstein action for gravity. (We note that one motivation for modifying Einstein's 
action is to attempt to do away with the construct of dark matter and/or dark energy. While it is 
probably too early to tell whether this can get rid of dark energy, it seems unlikely that this is a 
viable way of getting around the idea of dark matter.) 

For original papers written in languages other than English, we cite only an English translation, 
unless this does not exist. We only cite books that are in English. For books that have been 
reprinted we cite only the most recent printing of which we are aware. 

As a supplement to this review, we have compiled lists of key additional reference materials 
and links to Web resources that will be useful for those who want to learn more about this vast 
topic. These materials, available on our Web sit^j], include lists of more technical books (including 
standard textbooks on cosmology and related topics), historical and biographical references, less 
technical books and journal articles, and Web sites for major observatories and satellites. 



2. Foundations of the Big Bang Model 
2.1. General relativity and the expansion of the Universe 



Modern cosmology begins with lEinsteinl (|1917l ) where he applies his general relativity theory 
to cosmology. At this point in time our Galaxy, the Milky Way, was thought by most to be the Uni- 
verse. To make progress Einstein assumed the Universe was spa tially homog e neous an d isotr opic; 
this was enshrined as the "Copernican" cosmological principle bv lMilnd (|1933l ). iPeebled (|1993l . Sec. 
3) reviews the strong observational evidence for large-scale statistical isotropy; observational tests 
of homogeneity are not as straightforward. Einstein knew that the stars in the Milky Way moved 
rather slowly and decided, as everyone had done before him, that the Universe should not evolve in 
time. He could come up with a static solution of his equations if he introduced a new form of energy, 
now called the cosmqlogic al constant. It turns out that Einstein's static model is unstable. In the 
same year 
equations; 



de Sitter! (119171) fou nd the second cos mological solution of Einstein's general relativity 



Lemaitrd (119251 ) and iRobertsonl (119281 ) re-expressed this solution in the currentl y more 



familiar form of the exponentially expanding model used in the inflation picture. IWeyll (119231 ) noted 
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the importa nce of prescr i bing initial conditions such that the particle geodesies diverge from a point 
in the past. iFriedmannl (|1922l . [1924J), not bound by the desir e to have a stati c model, discovered 
the evolving homogeneous solu tions of Einstein 's equations; iLemaitrd (|1927l ) rediscovered these 
"Priedmann-Lemaitre" models. iRobertsonl (|1929i ) initiated the study of metric tensors of spatially 
homogeneous and isotropic spacetimes, and continuing study by him and A. G. Walker (in the mid 
1930's) led to the "Robertson- Walker" form of the metric tensor for homogeneous world models. 
Of course, in the evolving cosmological model solutions only observers at rest with respect to the 
expansion /contracti on see an iso tropic and homogen eous Universe; cosmology thus re-introduces 
preferred observers! iNorthl (|1990l ) and lLongairl (|2006l ) provide comprehensive historical reviews. 



See the standard cosmology textbooks for the modern formalism. 



2.2. Galaxy Redshift and Distance Measurements 



Meanwhile, with first success in 1912, Isiipheri (|l91# l finds that most of the "white spiral 
nebulae" (so-called because they have a continuum spectrum; what we now term spiral galaxies) 
emit light that is redshifted (we now know that the few, including M 31 (Andromeda) a nd some in the 
Virgo cluster, that emit blu e-shifted l i ght a re approaching us), and lEddingtonl (|1923l ) ide ntifies this 



with a redshift effect in the 



de Sitterl (|1917I ) model (not the cosmological redshift effect). iLemaitre 
(|1925l ) and IRobertsonl (|1928l ) derive Hubble's velocity-distance law v = Hqt (relating the galaxy's 



speed of recession v to its distance r from us, where Hq is the Hubble constant, the present value 
of the Hubble parameter) in the Friedmann-Lemaitre models. The velocity-distance Hubble law is 
a consequence of the cosmological principle, is exact, and implies that galaxies further aw ay than 
the cu rrent Hubble dist ance th = c/Hn are moving away faster than the speed of light c. iHubbk 



(192 



uses 



Leavitt's ( Leavitt 1912 : Johnson 2005 ^1 quantitative Cepheid variable star period- 



luminosity relation to establish th at M31 and M33 are far away (confirming the earlier somewhat 
tentative conclusion of I0piklll922l ). and does this for more galaxies, conclusively establishing that 
the white nebulae are other galaxies outside our Milky Way galaxy (there was some other earlier 
observational evidence for this position but Hubble's work is what convinces people). Hubble gets 
Humason (middle school dropout and one time muleskinner and janitor) to re-measure some Slipher 



2 Although the "canals" on Mars are not really canals, they had an indirect but profound influence on cosmology. 
Percival Lowell built Lowell Observatory to study the Solar System and Mars in particular, and closely directed the 
research of his staff. Slipher was instructed to study M31 and the other white nebulae under the hope that they were 
proto-solar-systems. 

3 Duncan had earlier found evidence for variable stars in M33, the spiral galaxy in Triangulum. 

4 Leavitt published a preliminary result in 1908 and Hertzsprung and Shapley helped develop the relation, but 
it would be another 4 decades (1952) before a reasonably accurate version became available (which led to a drastic 
revision of the distance scale). 
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spectra and measure more spectra, and iHubbld ( 1929l j 5 l establishes Hubble's redshift-distance law 
cz = Hqt, where the redshift z is the fractional change in the wavelength of the spectral line under 
study (although in the paper Hubble calls cz velocity and does not mention redshift). The redshift- 
distance Hubble law is an approxim ation to the velocity-distance law, valid only o n short distances 



and a t low redshif t s. iNorthl (|1990l ) provides a comprehensive historical review; iBerendzen et al 



(|1976l ) and ISmithl (jl982i ) are more ac cessible h istori cal summaries. See the standard cosmology 
textbooks for the modern formalism. iBranchl (|1998l ) discusses th e use of typ e la su pernovae as 



standard candles for meas uring the Hubb le constant. See Fig. 1 o f iLeibundeutl (120011) for a recent 



plot o f the Hubble law. iHarrisonl (|1993l ). iDavis &: Lineweaverl (|2004l ). and iLineweaver Davis 



(|2005l ) provide pedagogical discussions of issues related to galaxies moving away faster than the 
speed of light. 



2.3. The Hot Early Universe and Nucleosynthesis 



As one looks out further in space (and so back in time, because light travels at finite speed) 
wavelengths of electromagnetic radiation we receive now have been redshifted further by the expan- 
sion and so Wien's law tells us (from the blackbody CMB) that the temperature was higher in the 
past. The younger Universe was a hotter, denser place. Lemaitre ("the father of the Big Bang") 
emphasized the importance of accounting for the rest of known physics in the general relativistic 
cosmological models. 

Early work on explaining the astrophysically observed abundances of elements assumed that 
they were a consequence of rapid thermal equilibrium reactions and that a rapidly falling tem- 
perature froze the equilibrium abundances. Tolman, Suzuki, von Weizsacker, and others in the 
1920's and 1930's argued that the observed helium-hydrogen ratio in this scenario required that 
at some point the tempe r ature had been at least 10 9 K (and possibly as much as 10 11 K). 
Chandrasekhar Henrichl (119421 ) performed the first detailed, correct equilibrium computation 
and concluded that no single set of temperature and density values can accommodate all the ob- 
served abunda nces; they suggested that it would be useful to consider a non-equilibrium process. 
Gamowl (| 19461 ). building on his earlier work, makes the crucial point that in the Big Bang Model 
"the conditions necessary for rapid nuclear reactions were existing only for a very short time, so 
that it may be quite dangerous to speak about an equilibrium state", i.e., the Big Bang was the 
place to look for this non-equilibrium process. 



Gamowl (| 19481 ). a student of Friedmann, and lAlpherl (|1948l ) . a student of Gamow, estimated the 



radiation (photon) temperature at nucleosynthesis, and from the Stefan-Boltzmann law for black- 
body radiation noted that the energy budget of the Universe must then have been dominated by 



5 In the mid 1920's Lundmark and Stromberg had already noted that more distant galaxies seemed to have spectra 
that were more redshifted. 
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radiation. iGamowl (| 19481 ) evolved the radiation to the much later epoch of matter-radiation equal- 
ity (the matter and radiation energy densities evolve in different ways and this i s the time at which 
both had the same magnitude), a concept also introduced by Gamow, while lAlpher <fe Herman 
(| 19481 ) predicted a residual CMB radiation at the present time from nucleosynthesis and estimated 
its present temperature to be 5 K (because the zero-redshift baryon density was not reliably known 
then, it is s omewhat of a c oincidence that this temperature estimate is close to the observed mod- 
ern value). lHayashil (ll950l ) pointed out that at temperatures about 10 times higher than during 
nucleosynthesis rapid weak interactions lead to a thermal equilibrium abundance ratio of neutrons 
and protons determined by the neutron-proton mass difference, which becomes frozen in as the 
expansion decreases the temperature, thus establishing the initial conditions for nucleosynthesis. 
This is fortunate, in that an understanding of higher energy physics is not needed to make firm nu- 
cleosynthesis predictions; this is also unfortunate, because element abundance observations cannot 
be used to probe higher energy physics. 



Alpher et al.1 (|1953l ) conclude the early period of the standard model of nucleosynthesis. By 
this point it was clear that initial hopes to explain all observed abundances in this manner must fail, 
because of the lack of stable nuclei at mass numbers 5 and 8 and because as the temperature drops 
with the expansion it becomes more difficult to penetrate the Coulomb barriers. Cosmological 
nucleosynthesis can only generate the light elements and the heavier elements are generated from 
these light elements by further processing in the stars. 



Zel'dovichl ( 1963al ) and ISmirnovl ( 1964 ) noted that the 4 He and D abundances are sensitive 



to the baryon density : the observed abundances can be used to constrain the baryon density. 
Hovle fe Tavlerl ( 1964 ) carried out a detailed computation of the 4 He abundance and on comparing 
to measurements concluded "most, if not all, of the material of our .... Universe, has been 'cooked' 
to a temperature in excess of 10 10 K" . They were the first to note that the observed light element 
abundances were sensitive to the expansion rate during nucleosynthesis and that this could constrain 
new physics at that epoch (especially the number of light, relativistic, neutrino families). 



After Penzias and Wilson measured the CMB (see below), iPeeblesI (]1966al )bl) computed the 
abundances of D, 3 He, and 4 He, and their dependence on, among ot her things, the baryon density 
and the expansion rate during nucleosynthesis. The monumental IWagoner et al. I (jl967h paper 
established the ground r ules f or future work. For a history of these developments see pp. 125-128 
and 240-241 of lPeeblesI (|l97lL the articles by Al pher and Herman and Wago ner on pp. 129-157 



and 159-185, respectively, of lBertotti et al.1 (|1990h . and Ch. 3 and Sec. 7.2 of lKraghl (|1996l ). See 
the standard cosmology textbooks for the modern formalism. Accurate abundance predictions 
require involved numeri cal analysis; on the o t her hand pedagogy coul d benefit from approximate 



semi-analytical models ([Bernstein et al.lll989l : lEsmailzadeh et al.lll991 



In the simplest nucleosynthesis scenario, the baryon density estimated from the observed 
D abundance is consistent with that estimated from WMAP CMB anisotropy data, and higher 
than that estimated from the 4 He and 7 Li abundances. This is further discussed in Sec. 17.21 
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Field fc Sarkarl (|2006l ). ICvburtl (|2004l ) . and ISteigmanl (|200q ) are recent reviews of nucleosynthesis. 



In addition to residual CMB radiation, there is also a residual neutrino background. Above a 
temperature of about 10 10 K the CMB photons have enough energy to produce a thermal equilib- 
rium abundance of neutrinos. Below this temperature the neutrinos decouple and freely expand, 
resulting in about 300 neutrinos per cubic centimeter now (with three families, and this number 
also includes antineutrinos), at a temperature of about 2 K, lower than t h at of the CMB be cause 
electron-positron annihilation heats the CMB a little. See iDolgovi (|2002l ). iHannestadl (|2006l ). and 
the more recent textbooks cited below for more detailed discussions of the (as yet undetected) 
neutrino background. We touch on neutrinos again in Sec. 15.11 



2.4. Theory and Observations of the CMB 



The CMB radiation contributes of the order of 1% of the static or "snow" seen when switching 
between channels on a television with a conventional VHF antenna; it is therefore not surprising 
that it had been det ected a number o f times before its 1965 discovery/identification. For instance, 
it is now known that iMcKellari (Il94ll ) deduced a CMB temperature of 2.3 K at a wavelength of 2.6 
mm by estimating the ratio of populations in the first excited rotational and ground states of the 
interstellar cyanogen (CN) molecule (determined from absorption line measurements of Adams). It 
is now also known that the discrepancy of 3.3 K between the measured and expected temperature 
of the Bell L abs horn ant enna (for communicating with the Echo I satellite) at a wavelength of 12.5 
cm found by Ohm (1961) i s due to the CMB. Ohm also notes that an earlier measurement with this 
telescope (jDeGrasse et al.lll959i ) ascribes a temperature of 2 ± 1 K to back and side lobe pick up, 
that this is "...temperature not otherwise accounted for...", and that "it is somewhat larger than 
the calculated temperature expected" . Of course, McKellar had the misfortune of performing his 
analyses well before Gamow and collaborators had laid the nucleosynthesis foundations that would 
eventually explain the CN measurements (and allow the CMB interpretation) and Ohm properly 
did not overly stress the discrepancy beyond its weak statistical significance. 



While Alpher and Herman [e.g., pp. 114-115 of lAlpher fc Herman! (120011 ) and p. 130 of lWeinberg 



(|1993l )] privately raised the issue of searching for the CMB, and Hoyle c a me close to correct ly 
expla ining McKellar's CN measurements [se e pp. 345-346 oflKra eh (199601, IZel'dovi ch ( 196 3allbl. 
1965, p. 491, p. 89, and p. 315, respectively). iDoroshkevich &: Novikovi (|1964l ) and lDicke fc Peebles 
(jl965 . p. 448) are the first published discussions of possible observational consequences of the 
(then still hypothetical) CMB in the present Universe. The relevant discussions of Zel'dovich and 
Doroshkevich & Novikov are motivated by the same nucleosynthesis considerations that motivated 
Gamow and collaborators; Dicke and Peebles favored an oscillating Universe and needed a way to 
destroy hea yy elements from the previous cycle and so po s tulate d an ini tial hotter s tage in each 
cycle. Both IDoroshkevich &: Novikovi ([1964]) and[Zcl'dovich (1965) refer to lOhml (|196ll ) but neither 
appear to notice Ohm's 3.3 K discrepancy; in fact IZel'dovichl ()1965l ) (incorrectly) argues that Ohm 
constrains the temperature to be less than 1 K and given the observed helium abundance this rules 
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out the hot Big Bang Model! 



Working with the same antenna as Ohm, using the Dicke switching technique to compare the 
antenna temperature to a liquid heliu m load at a known temp erature, and paying very careful 
attention to possible systematic effects. iPenzias k. Wilson! ()1965l ) measure the excess temperature 
to be 3.5 i 1 K at 7.35 cm wavelength; Dicke et al. ( 19651 ) identify this as the CMB radiation left 
over from the hot Big Bang. 

The CMB is the dominant component of the radiation density of the Universe, with a density 
now of about 400 CMB photons per cubic centimeter at a temperature of about 2.7 K now. As noted 
in the previous subsection, observed light element abundances in conjunction with nucleosynthesis 
theory allows for constraints on the density of baryonic matter. Thus, there are a few billion CMB 
photons for every baryon; the CMB photons carry most of the cosmological entropy. 

To date there is n o observat i onal indication of any deviation of the CMB spectrum from a 
Planckian blackbody. iPartridgd (|1995l ) reviews early measurements of the CMB spectru m. A 
definitive observation of the CMB spectrum was made by COBE (see iGush et al.1 (Il990l ) for a 
contemporaneous ro cket-based measure ment), which measured a temperature of 2.725 ± 0.002 K 
(95 % confidence) (jMather et al.1 Il999l ) and 95 % confidence upper limits on possible spectral 
distortions: |/i| < 9 x 10~ 5 for the chemical potential of early (10 5 < z < 3 x 10 6 ) energy releas e 
and |y| < 1.5 x 10~ 5 for Comptonization of the spectrum at later times ( Fixsen et al.1 1 19961 ). 



Wright et al.l (|1994l ) shows that these constraints strongly rule out many alternatives to the Big 
Bang Model, including the steady state model and explosive galaxy formation. 



Anisotropy of the CMB temperature, first detected by COBE (ISmoot et al.l Il992l ). reveals 
important features of the formation and evolution of structure in the Universe. A small dipole 
anisotropy (discovered in the late 1960's and early 1970's by Conklin and Henry and confirmed by 
Corey and Wilkinson as well as Smoot, Gorenstein, and Muller) is caused by our peculiar motion; 
the CMB establishes a preferred reference frame. Higher multipole anisotropies in the CMB reflect 
the effect of primordial inhomogeneities on structure at the epoch of recombination and more 
complex astrophysical effects along the past light cone that alter this primordial anisotropy. We 
discuss these anisotropies, as well as the recently-detected polarization anisotropy of the CMB 
in Sec. 15.21 The anisotropy signal from the recombination epoch allows precise estimation of 
cosmological parameters (see Sec. 17.2ft . 

In addition to refere nces cited a bove, Dicke! ( 197C . pp. 64-70), Wilson ( 1983 ). Wilkinson Peebles 



(jl990l ). IPartridgd (|1995l . Ch. 2), and iKraghl (|1996l . Sec. 7.2) rev iew the history. For the modern for 
malism see the more recent standard cosmology textbooks and iKamionkowski Kosowskvl ([1999) 
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2.5. Challenges for the Big Bang Model 



Since the Universe is now expanding, at earlier times it was denser and hotter. A naive 
extrapolation leads to a (mathematical) singularity at the beginning, with infinite density and 
temperature, at the initial instant of time, and over all space. This naive extrapolation is unjustified 
since the model used to derive it breaks down physically before the mathematical singularity is 
reached. Deriving the correct equations of motion for the very early Universe is an important area 
of current research. While there has been much work, there is as yet no predictive model that 
unifies gravity and quantum mechanics — and this appears essential for an understanding of the 
very early Universe, because as one goes back in time the gravitational expansion of the Universe 
implies that large cosmological length scales now correspond to tiny quantum mechanical scales 
in the very very early Universe. There is a small but active group of workers who believe that 
only a resolution of this issue (i.e., the derivation of a full quantum theory of gravity) will allow 
for progress on the modeling of the very early Universe. But most others, perhaps inspired by the 
wonderful successes of particle physics models that have successfully described shorter and shorter 
distance physics, now believe that it is important to try to solve some of the "problems" of the Big 
Bang Model by attempting to model the cosmophysical world at an energy density higher than is 
probed by nucleosynthesis and other lower redshift physics, but still well below the Planck energy 
density where quantum gravitational effects are important. This is the approach we take in the 
following discussion, by focusing on "problems" that could be resolved below the Planck density. 
Whether Nature has chosen this path is as yet unclear, but at least the simplest versions of the 
inflation scenario (discussed in the next section) are compatible with current observations and will 
likely be well tested by data acquired within this decade. 

Assuming just nonrelativistic matter and radiation (CMB and neutrinos) in order of magnitude 
agreement with observations, the distance over which causal contact is possible grows with the age 
of the Universe. That is, if one assumes that in this model the cosmological principle is now valid 
because of "initial conditions" at an earlier time, then those initial conditions must be imposed 
over distances larger than the distance over which causal communication was possible. (And maybe 
this is what a quantum theory of gravitation will do for cosmology, but in the spirit of the earlier 
discussion we will view this as a "problem" of the Big Ba n g Mo del that should be resolved by 
physics at energies below the Planck scale.) lAlpher et al.l (|1953l . p. 1349) contains the earliest 
remarks (in passing) th at we are aware of about this particle horizon probl em. The termin ology 



is due to iRindlerl (119561 ) which is an early discussion o f horizons in g eneral. 



Harrison (1968) also 



mentions the particle horizon problem in passing, but iMcCreal (|1968l ) and iMisneri (|1969l ) contain 
the first clear statements we are aware of, with Misner stating "These Robertson- Walker models 
therefore give no insight into why the observed microwave radiation from widely differ ent angles in 
the sky has ... very precisely ... the sam e temperature" . Other early discu s sions are in iDicke ( 197C , 



Misner et al 



(1973 



p. 61), [D oroshke vich &: Novikovi (|l970l ). and the text books of IWeinberel (|1972l . pp. 525-526) and 



pp. 815-816), T his issue was dis c ussed in many papers and books starting in 



the early 1970's, but the celebrated IDicke Peebles! (|1979l ) review is often credited with drawing 
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prime-time attention to the particle horizon "problem" . 

The large entropy of the Universe (as discussed above, there are now a few billion CMB photons 
for every baryon) poses another puzzle. When the Universe was younger and hotter there had to 
have been a thermal distribution of particles and antiparticles and, as the Universe expanded and 
cooled, particles and antiparticles annihilated into photons, resulting in the current abundance of 
CMB photons and baryons. Given the lack of a significant amount of antibaryons now, and the 
large photon to baryon ratio now, at early times there must have been a slight (a part in a few 
billion) excess of baryons over antibaryons. We return to this issue in the next section. 



Inflation 



It is possible to trace a thread in the particle horizon problem tapestry back to the singularity 
issue and early discussions of Einstein, Lemaitre, and others who viewed the singularity as arising 
from the unjustified assumption of exact isotropy, and led to the intensive study of homogeneous 
but anisotropic cosmological models in the late 1960's and early 1970's. These attempts failed to 
tame the singularity but did draw attention to isotropy and the particle horizon problem of the 
standard Big Bang Model. It is interesting that this singularity issue also drove the development of 
the steady state picture, which in its earliest version was just a de Sitter model. While observations 
soon killed off the original steady state model (a more recent variant, the q uasi-steady state model, 



can be adjusted to accommodate the data, see, e.g., iNarlikar et al.l 12003 ) , the idea of a po ssible 



early, pre-Big-Bang, nonsingular de Sitter epoch thrived. It appears that iBrout et al.l (|1978l ) were 
the first to note that such a cosmological model was free of a particle horizon. However, they do 
not seem to make the connection that this could allow for isotropy by ensuring that points well 
separated now shared some c ommon ey ents in the past and thus causal physics could in principle 
make the Universe isotropic. IZed ()1980l ) noted that if one modifies the early Universe by speeding 
up the expansion rate enough over the expansion rate during the radiation dominated epoch, the 
particle horizon problem is resolved (but he does not go to the exponentially expanding de Sitter 
solution characteristic of the early inflation scenario) . 



Sato ( 


1981allbh. 


Kazanae ( 


1980), and 


Guth ( 


1981 



be crucial) point that during a phase transition at very high temperature in grand unified models 
it is possible for the grand unified Higgs scalar field energy density to behave like a cosmological 
constant, driving a de Sitter exponential cosmological expansion, which results in a particle- horizon- 
free cosmological model. And the tremendous expansion during the de Sitter epoch will smooth 
out wrinkles in the matter dist r ibutio n, by stretching them to very large scales, an effect alluded 
to earlier by iHovle &: Narlikarl ()1962j ) in the context of the steady state model, which could re- 
sult in an isotropic U nivers e now, provided the init i al wr inkles satisfy certain conditions. See 
Ellis Stoegerl (|1988l ) and INarlikar fc Padmanabhanl (|199ll ) for caveats and criticism. Of course, 
to get the inflationary expansion started requires a large enough, smooth enough initial patch. 
The contemporary explanation appeals to probability: loosely, such a patch will exist somewhere 
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and inflation will start there. In addition the initial conditions issue is not completely resolved 
by inflation, only greatly alleviated; since inflation stretches initially small length scales to length 
scales of contemporary cosmological interest, the cosmological principle requires that there not be 
very large irregularities on very small length scales in the very early Universe. This could be a clue 
to what might be needed from a model of very high energy, pre- inflation, physics. For reviews of 
inflation see the more recent standard cosmology and astroparticle physics textbooks. 



Building on ideas of Brout, Englert, and collaborators, iGottl (|1982l ) noted that it was possible 
to have inflation result in a cosmological model with open spatial hypersurfaces at the present time, 
in contrast to the Sato-Kazanas-Guth discussion that focused on flat spatial hypersurfaces. This 
open-bubble inflation model, in which the observable part of the contemporary Universe resides 
inside a bubble nucleated (because of a small upward "bump" in the potential energy density 
function) between two distinct epochs of inflation, is a clear counter-example to the oft-repeated 
(but incorrect) claim that infla t ion explain s why the Universe appears to have negligible space 
curvature. See iRatra &: Peebles! (|1994l . Il995l ) for a more detailed discussion of this model. 



The open-bubble inflation model was the first consistent inflation model. Unfortunately for 
the Guth model, as the phase transition completes and one hopes to have a smooth transition to 
the more familiar radiation-dominated expansion of the hot Big Bang Model, one finds that the 
potential in the Guth model results in many small bubbles forming with most of energy density 
residing in the bubble walls. In this model the Universe at the end of inflation was very inhomo- 
geneous because the bubble collisions were not r apid enough t o the rmalize the bubble wall energy 
density (i.e., the bubbles did not "percolate"). iLindd (|1982l ) and lAlbrecht Sz Steinhardtl (|1982l ) 
used a specific potential energy density function for the Higgs field in a grand unified model and 
implemented Gott's scenario in the Sato-Kazanas-Guth picture, except they argued that the sec- 
ond epoch of inflation lasts much longer than Gott envisaged and so stretches the bubble to length 
scales much larger than the currently observable part of the Universe, thus resulting in flat spatial 
hypersurfaces now. The great advantage of the Gott scenario is that it uses the first epoch of infla- 
tion to resolve the particle horizon/homogeneity p roblem and so this prob lem does not constrain 
the am ount of inflation after the bubble nucleates. iBrout et al.l (|1978l ) and lColeman &: De Luccia 
(|1980l ) note that symmetry forces the nucleating bubble to have an open geometry, and this is why 
inflation requires open spatial hypersurfaces, but with significant inflation after bubble nucleation 
the radius of curvature of these hypersurfaces can be huge. Thus the amount of space curvature in 
the contemporary Universe is a function of the amount of inflation after bubble nucleation, and it 
is now widely accepted that observational data (as discussed below in Sec. I7.2p are consistent with 
an insignificant amount of space curvature and thus significant inflation after bubble nucleation. 

It is well known that phase transitions can create topological defects. Grand unified phase 
transitions are no exception and often create monopoles and other topological defects. If the 
Universe is also inflating through this phase transition then the density of such topological defects 
can be reduced to levels consistent with the observations. This is not another argument in support 
of inflation, although it is often claimed to be: it is just a way of using inflation to make viable a 
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grand unified theory that is otherwise observationally inconsistent. 



One major motivation for grand unification is that it allows for a possible explanation of the 
observed excess of matter over antimatter (or the baryon excess) mentioned in the previous section. 
There are other possible explanations of how this baryon excess might have come about. One much 
discussed alternative is the possibility of forming it at the much lower temperature electroweak phase 
transition, through a non-perturbative process, but this might raise particle horizon or homogeneity 
issues. However, at prese nt there is no convincing, numerically satisfying exp l anation of the b aryon 
excess , from any process. lOuinn Witherelll (|1998l ). iDine fc Kusenkd (|2004l ) . iTroddenl (|2004l ). and 
Clind (|2006l ) review models now under discussion for generating the baryon excess. 



At the end of inflation, as the phase transition completes and the Universe is said to reheat, 
one expects the generation of matter and radiation as the Universe makes the transition from rapid 
inflationary expansion to the more sedate radiation-dominated expansion of the hot Big Bang 
Model. This is an area of ongoing research and it would be useful to have a convincing, numerically 
satisfying model of this epoch. The baryon excess might be generated during this reheating process. 

While great effort has been devoted to inflation, resulting in a huge number of different models, 
at the present stage of development inflation is a very interesting general scenario desperately in 
need of a more precise and more convincing very high energy particle physics based realization. 
As far as large-scale cosmology is concerned, inflation in its simplest form is modeled by a scalar 
field (the inflaton) whose potential energy density satisfies certain properties that result in a rapid 
enough cosmological expansion at early times. It is interesting that cosmological observations 
within this decade might firm up this model of the very early Universe based on very high energy 
physics before particle physicists do so. For reviews see the more recent standard cosmology and 
astroparticle physics textbooks. 

Assuming an early epoch of inflation, the cumulative effect of the expansion of the Universe 
from then to the present means that contemporary cosmological length scales (e.g., the length 
scale that characterizes the present galaxy distribution) correspond to very tiny length scales dur- 
ing inflation, so tiny that quantum-mechanical zero-point fluctuations must be considered in any 
discussion involving physics on these length scales. 

As mentioned above, the idea of an early de-Sitter-like expansion epoch, pre-Big-Bang, was 
discussed in the 1970's, as a possible way of taming the initial singularity. While this de Sitter 
epoch was typically placed at very high energy, it differs significantly from the inflation scenario 
in that it was not driven by a scalar field potential energy density. Nevertheless because it was at 
energies close to the Planck energy there were many discussions of quantum mechanical fluctuations 
in de Sitter spacetime in the 1970's. 

In the inflation case quantum mechanics i ntroduces addit i onal fluctuations, th e zer o-point fluc- 
tuations in the scalar field. This was noted by Hawking) (|1982l ). IStarobinskvl (|1982l ). and lGuth fc Pil 
(|l982h . and further studied bv bardeen et all JlQSsK . T or a discussion o f scalar field quantum 
fluctuations in de Sitter spacetime and their consequences see iRatral (| 19851 ) . iFischler et al.l (|1985l ) 



-13- 



use the Dirac-Wheeler-DeWitt formalism to consistently semi-classically quantize both gravitation 
and the scalar field about a de Sitter background, and carry through a computation of the power 
spectrum of zero-point fluctuations. The simplest inflation models have a weakly coupled scalar 
field and so a linear perturbation theory computation suffices. The fluctuations obey Gaussian 
statistics and so can be completely characterized by their two-point correlation function or equiv- 
alently their po wer spectrum. Inflat ion models that give non-Gaussian fluctuations are possible 
(for a review see lBartolo et al.l (|2004l )). but the observations do not yet demand this, being almost 
completely consistent with Gaussianity (see discussion in Sec. l5.2l below). The simplest models give 
adiabatic or curvature (scalar) fluctuations; these are what result from adiabatically compressing 
or decompressing parts of an exactly spatially homogeneous Universe. More complicated models 
of inflation can produ ce fluctuations that break adi abaticity, such as (tenso r ) gravitationa l waves 



(jRubakov et al.lll982l ) and (vector) magnetic fields ([Turner k, Widrow 



1988 



Ratralll992al ). which 



might have interesting observational consequences (see Sees. [5TT1 and [ST21 below). 

The power spectrum of energy density fluctuations depends on the model for inflation. If the 
scalar field potential energy density during inflation is close to flat and dominates the scalar field 
energy density, the scale factor grows exponentially with time (this is the de Sitter model), and 
after inflation but at high redshift the power spectrum of (scalar) mass-energy density fluctuations 
with wavenumber magnitude k is proportional to k, or scale invariant, on all interesting length 
scales, i.e., curv ature fluctuations diverge only as logfc. This was noted in the early 1980's for the 



inflation model (|Hawkinglll982l ; IStarobinskv 



1982 



Guth fc Pill 19821 ) , although the virtues of a scale - 



invariant spec t rum were emphasized in the early 1970's, well before inflation, by iHarrisonl (|l970l ). 
Peebles Yul (|l970l ). and lZel'dovichl (|l972l ). When the scalar field potential energy density is such 
that the scalar field kinetic energy density is also significant during inflation a more general spectrum 
proportional to k n can result (where the spectral index n depends on the slope of the potential 
energy density during inflation); for n ^ 1 the spectrum is said to be tilted (lAbbott fc Wisdll984l : 



Lucchin fc Matarrese 



19851 ; lRatralll992bl ). Current observations appear to be reasonably well fit by 



n 



1. More complicated, non-power-law spectra are also possible. 



We continue this discussion of fluctuations in Sec. [5] below. 



4. Dark Matter and Dark Energy 

Most cosmologists are of the firm opinion that observations indicate the energy budget of the 
contemporary Universe is dominated by dark energy, with the next most significant contributor 
being dark matter, and with ordinary baryonic matter in a distant third place. Dark energy and 
dark matter are hypothetical constructs generated to explain observational data, and the current 
model provides a good, but not perfect, explanation of contemporary cosmological observations. 
However, dark energy and dark matter have not been directly detected (in the lab or elsewhere). 
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Hubbld (1926J_| presented the first systematic estimate of masses of the luminous part of 
galaxies (based on studying the motion of stars in galaxies), as well as an estimate of the mass 
density of the Universe (using counts of galaxies in conjunction with the estimated masses of 
galaxies) . 

Under similar assumptions (the validity of Newton's second law of motion and Newton's 
inverse-square law of gravitation, and that the large-scale structure under investigation is in gravita- 
tional equilibrium) . IZwickvl (j 19331 ). in perhaps one of the most significant discoveries of the previous 
century, found that galaxies in the Coma cluster of galaxies were moving with surprisingly high 
speeds. In modern terms, this indicates a Coma cluster mass density at least an order of magnitude 
greater than what would be expected from spreading the mass associated with the luminous parts 
of the galaxies in the Coma cluster over the volume of the cluster. Zwicky's measurements probe 
larger length scales than Hubble's and so might be detecting mass that lies outside the luminous 
parts of the galaxies, i.e., mass that does not shine, or dark matter. Ordinary baryonic matter is 
largely nonrelativistic in the contemporary Universe and hence would be pulled in by the gravita- 
tional field of the cluster. Nucleosynthesis and CMB anisotropy measurements constrain the mass 
density of ordinary baryonic matter, and modern data indicate that not only is the amount of grav- 
itating mass density detected in Zwicky-like observations significantly greater than what is shining, 
it is likely a factor of 3 to 5 times the mass density of ordinary baryonic matter. (It is also known 
that a large fraction of the expected baryonic matter can not si gnifican t ly shi ne.) ISmithl (| 19361 ) 
confirmed Zwicky's result, using Virgo cluster measurements, and Zwickv (1937)0 soon followed up 
with a more detailed paper. 



Babcock's Ph.D. thesis (|Babcocklll939l ) was the next major (in hindsight) development in the 
dark matter story. He measured the rotation speed of luminous objects in or near the disk of the 
Andromeda (M31) galaxy, out to a distance of almost 20 kpc from the center and found that the 
rotation speed was still rising, not exhibiting the 1/y/r Keplerian fall off with distance r from the 
center that would be expected if the mass distribution in M31 followed the distribution of the light. 
That is, Babcock found that the out er part of th e luminous part of M31 was dominated by matter 
that did not shine. Soon t hereafter Port ( 19401 not ed a similar result for the galaxy NGC 3115. 
Almost two decades later, van de Hulst et al.l (jl957 ) confirmed Babcock's result by using 21 cm 
wavelength observations of hydrogen gas clouds that extend beyond the luminous part of M31, 
finding a roughly flat rotation curve at the edge (no longer rising with distance as Babcock had 
found). While there was some earl y theoretical di s cussio n of this issue, the much more detailed 
M31 flat rotation curve measured bv lRubin Fordl (|1970h (Rubin was a student of Gamow) forced 



6 In this paper, among other things, Hubble also developed his galaxy classification scheme (of ellipticals, normal 
and barred spirals, and irregulars) and showed that the averaged large-scale galaxy distribution is spatially isotropic. 



Opik (1922) had earlier estimated the mass of M31. 



7 In this paper Zwicky also proposes the remarkable idea of using gravitational lensing of background objects by 
foreground clusters of galaxies to estimate cluster masses. 



-15- 



this dark matter into the limelight. 



Oth er early ind ications of dark matter came from measurements of the velocities of binar y 
galaxies (Pagd 1952) and the dynamics of our Local Group of galaxies (jKahn <fe Woltieii 11959). 



de Vaucouleurd (119691 ) and lArp Bertolal (|1969l) found that t he ell iptical galaxy M87 in the Virgo 



cluster had a faint mass-containing halo. lOstriker Peebles! ()1973l ) noted that one way of making 
the disk of a spi ral galaxy stab l e aga inst a bar-like instab il ity is to embed it in a massive halo, and 
soon thereafter lEinasto et al.l ()1974l ) and lOstriker et al.l ()1974l ) showed that this suggestion was 
consistent with the observational evidence. These early results have been confirmed by a number 
of different techniques, including measuring the X-ray temperature of hot gas in galaxy clusters 
(which is a probe of the gravitational potential — and the mass which generates it — felt by the 
gas), and measurements of gravitational lensing of background sources by galaxy clusters. See Sec. 
17.21 for further discussion of this. 



For reviews of dark matter see Sec. IV of iPeebled (|197ll ) (note the fascinating comment on p 



64 on the issue of dark matter in clusters: "This quantity" M/L or the mass to luminosity ratio "is 



(|l992h . lPeeblesl dl99.4 Sec. 18), andlanastcj (fcoosh . 



(1979). 


Trimble 


(1987 


), 


Ashman 



Much as van Ma anen's measu rements of large (but erroneous) rotation velocities for a number 
of galaxies prompted iJeansI (|1923l ) to consider a modification of Newton's inverse-square law for 
gravity such that the gravitational force fell off slower with distance on l arge d i stanc es, the large 
(but not erroneous) velocities measured by Zwicky and others prompted iFinzil (|1963l ) , and many 
since then, to consider modifications of the law of gravity. The current observational indicatio ns are 
that this is not a very viable alternative to the dark matter hypothesis ([Peebles Ratrall2003l . Sees. 
IV. A. 1 and IV.B.13). In some cases, modern high energy physics suggests possible motivations for 
modifications of the inverse square law on various length scales; this is beyond the scope of our 



review. 



Milgroml (|1983l . |2002j) proposes a related but alternate hypothesis: Newton's second law of 
motion is modified at low accelerations. This hypothesis - dubbed modified Newtonian dynamics 
(MOND) - does a remarkable job of fitting the flat rotation curves of spiral galaxies, but most 
who have cared to venture an informed opinion believe that it cannot do away completely with 
dark matter, especially in low-surface-brightness dwarf galaxies and rich clusters of galaxies. More 
importantly, the lack of a well motivated extension of the small-length-scale phenomenological 
MOND hypothesis that is applicable on large cosmological l ength scales greatl y hinders testing the 
hypothesis. For a recent attempt at suc h an extension see iBekensteinl (|2004j). For a preliminary 
sketch of cosmology i n thi s context see biaz-Rivera et al.l feood ). For a review of MOND see 
Sanders & McGaughl feood ). 



Most cosmologists are convinced that dark matter exists. Nucleosynthesis constraints indicate 
that most of the dark matter is not baryonic. (Not all baryons shine; for a review of options for dark 
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baryons see lCarrlll994l ). Galaxies are in general older than larger-scale structures (such as clusters); 
this indicates that the dark matter primeval velocity dispersion is small (for if it were large gravity 
would be able to overcome the corresponding pressure only on large mass — and so length — scales, 
first forming large-scale objects that fragment later into younger smaller-scale galaxies). Dark mat- 
ter with low primeval velocity dispersion is known as cold dark matter (CDM). More precisely, the 
CDM model assumes that most of the nonrelativistic matter-energy of the contemporary Universe 
is in the form of a gas of massive, non-baryonic, weakly-interacting particles with low p rimeval 
veloci t y dispersion. One rea son t hey must be weakly interacting is so they do not shine. iMuhoz 



(|2004l ) , iBertone et al.l (|2005l ) , and iBaltz (l2004ri revi e w par ticle physics dark matte r cand idates and 



prospects for experimental detection. iBond et al.l (|1982l ) and iBlumenthal et al.l (|1982l ) note the 
advantages of CDM and that modern high e nergy physics models provide plausible hypothetical 
candidates for these particles. iPeeblesI (|1982l ) casts the cosmological skeleton of the CDM model, 
emphasizing that in this model structure forms from the gravitational growth of primordial depar- 
tures from homogeneity that are Gaussian, a diabatic, and scale invar iant, consistent with what is 
expected from the simples t infla tion models. Blumentha l et al.1 (| 19841 ) is a first fleshing out of the 
CDM model. See IPeeblesI l|l99ah and Middle fc Lvthl ||200gI ) for textbook discussions of the CDM 
model. More details about this model, including possible problems, are given in Sec. [5] below. 



To set the numerical scale for cosmological mass densities, following Einstein &; de Sitter 



(jl932l ). one notes that the simplest Friedmann-Lemaitre model relevant to the contemporary Uni- 
verse is one with vanishing space curvature and with energy budget dominated by non-relativistic 
matter (and no cosmological constant). In this critical or Einstein-de Sitter case the Friedmann 
equation fixes the energy density of nonrelativistic matter for a given value of the Hubble con- 
stant. Cosmologists then define the mass-energy density parameter Q for each type of mass-energy 
(including that of the curvature of spatial hypersurfaces Ok, the cosmological constant and 
nonrelativistic matter Qm) as the ratio of that mass-energy density to the critical or Einstein-de 
Sitter model mass-energy density. The Friedmann equation implies that the mass-energy density 
parameters sum to unity. (In general the fi's are time dependent; in what follows numerical values 
for these parameters refer to the current epoch.) 

As discussed in Sec. 17.21 below, it has long been known that nonrelativistic matter (baryons 
and CDM) contributes about 25 or 30 % to the critical mass-energy density. After the development 
of the inflation picture for the very early Universe in the 1980's there was a wide-spread belief 
that space curvature could not contribute to the mass-energy budget (this is not necessary, as 
discussed above), and for this and a few other reasons (among others, the time scale problem 
arisin g from the large measured values of the Hubble constant and age of the Universe), iPeebles 
(|1984l ) proposed that Einstein's cosmological constant contributed the remaining 70 or 75 % of 
the mass-energy of the Universe. This picture was soon generalized to allow the possibility of a 
scalar- field energy density that is slowly varying in time and close to ho mogeneous in space, what 



is now called dark energy (jPeebles &: Ratra 



1988 



Ratra Peebles! Il988l ). As discussed in Sec. 17.21 



below, these models predict that the expansion of the Universe is now accelerating and, indeed, it 
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appears that this acceleration has been d etected at about the magnitude predicted in these models 



(jRiess et al.lll998l ; iPerlmutter et al.l ll999). Consistent with this, CMB anisotropy observations are 



consistent with flat spatial hypersurfaces, which in conjunction with the low mass-energy density 



are discussed in more detail in Sec. 17.21 below and in reviews ( Peebles fc Ratra 



2003 
2006 



Carroll 



2004; 



Sahni fe Starobinsky 



200 



Padmana bhanll2005l : IPerivolaropoulosI 120061 : 



(Peebles & Ratra 


2003: 


Steinhardt 


Copeland et al. 


2006 


; i 


N"obbenhuis 



The following sections flesh out this "standard model" of cosmology, elaborating on the model 
as well as describing the measurements and observations on which it is based. 



5. Growth of Structure 
5.1. Gravitational Instability and Microphysics in the Expanding Universe 

5.1.1. Gravitational Instability Theory from Newton Onwards 

The primary driver for the formation of large-scale structure in the Universe is gravitational 
instability. The detailed growth of structure depends on the nature of the initial fluctuations, the 
background cosmology, and the constituents of the mass-energy density, as causal physics influences 
the rate at which structure may grow on different scales. 

Newton, prompted by questions posed to him by Bentley, realized that a gas of randomly 
positioned massive particles interacting gravitationally in flat spacetime is unstable, and that as 
time progresses the mass density distribution grows increasingly more anisotropic and inhomoge- 
neous. Awareness of this instability led Newton to abandon his preference for a finite and bounded 
Universe of stars for one that is infinite and homogeneous on average (see discussion in Irlarrisonl 



2001 



); this was an early discussion of the cosmological principle. 



Jean si (| 19021 ) studied the stability of a spherical distribution of gravitating gas particles in flat 
spacetime, motivated by possible relevance to the process of star formation. He discovered that gas 
pressure prevents gravitational collapse on small spatial scales and gives rise to acoustic oscillations 
in the mass density inhomogeneity, as the pressure gradient and gravitational forces compete. On 
large scales the gravitational force dominates and mass density inhomogeneities grow exponentially 
with time. The length scale on which the two forces balance has come to be known as the Jeans 
length or the acoustic Hubble length c s /Ho, where c s is the speed of sound. 

On scales smaller than the Jeans length, adiabatic energy density perturbations oscillate as 
acoustic waves. On scales well below the Jeans length dissipative fluid effects (e.g., viscosity and 
radiation diffusion) must be accounted for. These effects remove energy from the acoustic waves, 
thus damping them. In an expanding Universe, damping is effective when the dissipation time scale 
is shorter than the expansion time scale, and the smallest length scale for which this is the case is 
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called the damping length. This is discussed in more detail below. 



5.1.2. Structure Growth in an Expanding Universe 



Study of gravitational instability in an evolving spacetime, appropriate for the expanding 
Universe, began with Lemaitre in the early 1930's. He pioneered two approaches, both of which are 
still in use: a "nonperturbative" approach based on a spherically symmetric solution of the Einstein 
equations (further developed by Dingle, Tolman, Bondi, and others and discussed in Sec. l5.3l below); 
and a "perturbative" approach in which one studies small departures from spatial homogeneity and 
isotropy evolving in homogeneous and isotropic background spacetimes. 

At early times, and up to the present epoch on sufficiently large scales, the growth of structure 
by gravitational instability is accurately described by linear perturbation theory. The growth of 
small density and velocity perturbations must take into account the effects of the expansion of the 
Universe. A fully relativistic theory must be employed to describe the growth of structure, because 
it is necessary to also describe the evolution of modes with wavelength larger than the Hubble 
length. In contrast, a Newtonian approximation is valid and used on smaller length scales. 



L if shit a (| 19461 ) laid the foundations of the general-relativistic perturbative approach to struc- 



ture formation. He linearized the Einstein and stress-energy conservation equations about a spa- 
tially homogeneous and isotropic Robertson- Walker background spacetime metric and decomposed 
the departures from homogeneity and isotropy into independently evolving spatial harmonics (the 
so-called scalar, vector, and tensor modes). Lifshitz treated matter as a fluid which is a good 
approximation when the underlying particle mean free path is small. He discovered that the vector 
transverse peculiar velocity (the peculiar velocity is the velocity that remains after subtracting 
off that due to the Hubble expansion) perturbation decays with time as a consequence of angu- 
lar momentum conservation and that the contemporary Universe could contain a residual tensor 
gravitational wave background left over from earlier times. 

Unlike the exponentially growing energy density irregularity that Jeans found in flat spacetime 
on large scales, Lifshitz found only a much slower power-law temporal growth, leading him to the 
incorrect conclusion that "gravitational instability is not the source of conden sation of matt er into 
separate nebulae". It was almost two decades before iNovikovl (| 19641 ) (but see lBonnorlll9571 . for an 
earlier hint) corrected this misunderstanding, noting that even with power-law growth there was 
more than enough time for inhomogeneities to grow, since they could do so even while they were 
on scales larger than the Hubble length rjj = c/Hq at early time. 

The approach to the theory of linear perturbations initiated by Lifshitz is based on a specific 
choice of spacetime coordina tes called synch ronous coo rdinates. This approach is dis c ussed in deta il 
in Sec. V (also see Sec. II) of lPeeblesI jl98ch . Sec. Ill of lZel'dovich fc Novikovi ( 1983h .iRatrJ (jl988h . 
and other standard cosmology and astroparticle physics textbooks. iBardeenl (|1980l ) (building on 
earlier work) recast the Lifshitz analysis in a coordinate-independent form, and this approach has 
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also become popular. For reviews of this approach see iMukhanov et al.l (|1992), as well as the 
standard textbooks. 

A useful for malism for l inear growth of density and velocity fields is given by the "ZeFdovich 
approximation" (jZel'dovichl I197Q : IShandarin &: Zel'dovichl Il989i : ISahni fc Colesl Il995l ) , based on 
anisotropic collapse and so "pancake" formation (a concept earlier discussed in the context of the 
initial singularity). This method accurately describes structure formation up to the epoch when 
nonlinearities become significant. Numerical simulations (see Sec. 15.41 below) of fully non-linear 
structure growth often employ the Zel'dovich approximation for setting the initial conditions of 
density and velocity. 



5.1.3. Space Curvature 



The evolution of the backgroun d spacetime influe n ces th e rate of growth of structure. An early 
example of this effect is seen in the iGamow fc Teller! (jl939l ) approximate generalization of Jeans' 
analysis to the expanding Universe, in particular to a model with open spatial hyper surf aces. At 
late times the dominant form of energy density in such a model is that due to the curvature of 
spatial hyper surf aces, because this redshifts away slower than the energy density in nonrelativistic 
matter. The gravitational instability growth rate is determined by the matter energy density, 
but the expansion rate becomes dominated by the space curvature. As a result, the Universe 
expands too f ast for inhomog eneities to grow and large-scale structure formation ceases. [A quarter 
century later, iPeebled (|1965l ) noted the importance of this effect.] This was the first example of an 
important and general phenomenon: a dominant spatially-homogeneous contributor to the energy 
density budget will prevent the growth of irregularity in matter. 



5.1.4- Dark Energy 



Matter perturbations also cannot grow when a cosmological constant or nearly homogeneous 
dark energy dominates. There is strong evidence that dark energy — perhaps in the form of 
Einstein's cosmological constant — currently contributes ~ 70 % of the mass-energy density of the 
universe. This dark energy was su b-dominant until recently, when it started slowing the rate of 
growth of structure ([Peebles! 1 19841 ). thus its effect on dynamical evolution is milder than that of 
space curvature. 



5.1.5. Radiation and its Interaction with Baryonic Matter 



Guyot Zel'dovichl (jl970l ) showed that a dominant homogeneous radiation background makes 
the Universe expand too fast to allow matter irregularities to start growing until the model becomes 
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matter dominated (when the radiation redshifts away). Because of this effect, as discussed next, 
the acoustic Hubble length at the epoch when the densities of matter and radiation are equal is 
an important scale for structure formation in the expanding Universe. This imprints a feature in 
the power spectrum of matter fluctuations on the scale of the acoustic Hubble length at matter- 
radiation equality that can be used to measure the cosmic density of non-relativistic matter. We 
return to this in Sec. 17.2} a related CMB anisotropy effect is discussed in the next subsection I5.2L 



Gamowl (j 19481 ) noted that at early times in the Big Bang Model radiation (which has large 



relativistic pressure) dominates over baryonic matter. In addition, at high temperature radiation 
and baryonic matter are strongly coupled by Thomson-Compton scattering and so behave like a 
single fluid. As a result of the large radiation pressure during this early epoch the Jeans or acoustic 
Hubble length is large and so gravitation al growth of inhomo geneity occurs only on large scales, 
with acoustic oscillations on small scales. iPeebles fc Yul ()1970l ) develop this picture. 

As the Universe cooled down below a temperature T ~ 3000 K at a redshift z ~ 10 3 , the 
radiation and baryons decoupled. Below this temperature proton nuclei can capture and retain free 
electrons to form electrically neutral hydrogen atoms — this process is called "recombination" - 
because fewer photons remained in the h i gh en ergy tail of the distribution with enough energy to 
disassociate the hydrogen atoms. IPeebles! (119681 ) and lZel'dovich et al.1 ( 196a ) perform an analysis of 
cosmological recombination, finding that at the "end" of recombination there were enough charged 
particles left over for the Universe to remain a good conductor all the way to the present. The 
finite time required for recombination results in a surface of non-zero thickness within which the 
decoupling of now-neutral baryons and photons occurs. The mean-free path for photons quickly 
grew, allowing the photons to travel (almost) freely, thus this "last-scattering surface" is the "initial" 
source of the observed CMB; it is an electromagnetically opaque "cosmic photosphere". See the 
standard cosmology textbooks for discussions of recombination. 

Decoupling leads to a fairly stee p drop in the p ressure of the baryon gas, and so a fairly steep 
decrease in the ba ryon Jeans length. IPeebles! (|1965l ) was developi ng this picture as the CMB was 
being discovered. IPeebles &: Dickd (119681 ) (also see IPeebles! 1 1 96 71 ) noted that the baryonic Jeans 
mass after decoupling is of the order of the mass of a typical globular cluster and so proposed that 
proto-globular-clusters were the first objects to gravitationally condense out of the primordial gas. 
This model would seem to predict the existence of extra-galactic globular clusters, objects that 
have not yet been observationally recognized. There are, however, dwarf galaxies of almost equally 
low mass, and we now also know that some globular clu sters are young and so globular clusters 
might form in more than one way (for a recent review see lBrodie &: Straderll2006l ). 



On scales smaller than the Jeans mass, dissipative effects become important and the ideal fluid 
approximation for radiation and baryonic matter is no longer accurate. As the Universe cools down 
towards recombination and decoupling, the photon mean free path grows and so photons diffuse out 
of more dense regions to less dense regions. As they diffuse the photons drag some of the baryons 
with them and so damp small-scale inhomogeneities in the photon-baryon fluid. This collisional 



-21 - 



damping — a consequence of Thomson-Co mpton s c atteri ng — is known as Sil k damping in the 



cosmological context; it was first studied by iMichid (|1969l 1 8 l. Peebles! (|1967l ). and ISilkl ()1968l ). The 



Silk damping scale is roughly that of a cluster of galaxies. 



5.1.6. Possible Matter Constituents 

If baryons were the only form of non-relativistic matter the density of matter would be so 
low that the Universe would remain radiation dominated until after recombination. The expansion 
rate would be too large for gravitational instability to cause inhomogeneity growth until matter 
starts to dominate well after last scattering. The short time allowed for the gravitational growth of 
inhomogeneity from the start of matter domination to today would require a large initial fluctuation 
amplitude to produce the observed large-scale structure. This scenario is ruled out by measurements 
of the anisotropy of the CMB which indicate that fluctuations in the baryons at decoupling are 
too small to have grown by gravitational instability into the structures seen today in the galaxy 
distribution. 

A solution to this puzzle is provided by dark matter, of the same type and quantity needed 
to explain gravitational interactions on galactic and cluster scales. Including this component of 
matter the Universe becomes matter dominated at a redshift comparable to, or even larger than, 
the redshift of last scattering. Because CDM does not directly couple to radiation, inhomogeneities 
in the distribution of CDM begin to grow as soon as the Universe becomes matter dominated. 
Growth in structure in the baryons, on scales small compared to the Hubble length, remains 
suppressed by Thomson-Compton scattering until recombination, after which baryons begin to 
gravitate toward the potential wells of dark matter and the baryon fluctuation amplitude quickly 
grows. Thus, the low observed CMB anisotropy is reconciled with observed large-scale structure. 
(The CMB, while not directly coupled to the CDM, feels the gravitational potential fluctuations of 
the CDM. Consequently, measurements of the CMB anisotropy probe the CDM distribution.) This 
is an independent, although model-dependent and indirect, argument for the existence of CDM. 

As mentioned above in Sec. 12.3} the Universe also contains low mass neutrinos (precise masses 
are not yet known). These neutrinos are relativistic and weakly coupled (nearly collisionless) and 
so have a very long mean free path or free-streaming length. Consequently, they must be described 
by a distribution function, not a fluid. Because they are relativistic they have a large Jeans mass 
and gravitational instability is effective at collecting them only on very large scales, i.e., low mass 
neutrinos suppress power on small and intermediate length scales. This effect makes it possible 



to observationally probe th ese particles with cosmological measurements (jElgarOv fc Lahavil2005l ; 



Lesgourgues fc Pastor! 120061 ) 



This is a version of a manuscript submitted to the Astrophysical Journal on September 1, 1967, and only 
minimally revised (in response to the referee's suggestions) before the author died. 
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5.1.7. Free Streaming 

Thus, the properties of dark matter are reflected in the spectrum of densi ty fluctuations be cause 



scales smaller than the free-streaming scale of massive particles are damped (jBond et al.lll980l ) . For 
hot dark matter (HDM), e.g., neutrinos, the free-streaming scale is larger than the Hubble length 
at matter-radiation equality, hence the spectrum retains only large-scale power. In such a "top- 
down" scenario, superclusters form first, then fragment into smaller structures including clusters of 
galaxies and individual galaxies, as first discussed by Zel'dovich and collaborators. The top-down 
model was inspired by experimental suggestions (now known to be incorrect) that massive neutrinos 
could comprise the nonbaryonic dark matter, and by an early (also now known to be incorrect) 
interpretation of observational data on superclusters and voids (see Sec. !6.2l below) that postulated 
that these were the basic organizational blocks for large-scale structure. It predicts that smaller 
scale structure (e.g., galaxies) is younger than larger scale structure (e.g., superclusters), contrary 
to current observational indications. In fact, these observational constraints on the evolution of 
structure constrain the a mount of HDM neutrino matter-energy density and so neutrino masses 



(jKahniashvili et al.1 120051 ) . Cosmological observations provide the best (model-dependent) upper 



limits on neutrino masses. 

For CDM, e.g., weakly interacting massive particles (WIMPs), the free-streaming scale is 
negligible for cosmological purposes. This "bottom-up" or "hierarchical" scenario, pioneered by 
Peebles and collaborators, begins with the formation of bound objects on small scales that aggregate 
into larger structures, thus galaxies result from mergers of sub-galaxies, with superclusters being 
the latest structures to form. This is in better agreement with the observational data. See Sec. H] 
for more details on this model. 



5.1.8. Initial Density Perturbations and the Transfer Function 



The current standard model for structure formation assumes that structure in the Universe 
arose primarily from gravitational amplification of infinitesimal scalar density perturbations in 
the early Univers e. The process es listed in this section modify these initial in homogeneities. Re - 
~~3( 



views are given m 



Peebles! (I1980I. Sec. VUZel'dovich k Novikovi (11981 Sec. IIP . lEfstathioul (1 1990') 



Kolb Turned (|1990l . Ch. 9) . IPadmanabhanl (1 1993 . Ch. 4 UDekel fc Ostrikerl (|1999L and 
|2005|, Part U). 



Mukhanov 



As discussed in Sees. [H and 15.21 observations to date are consistent with primordial fluctua- 
tions that are Gaussian random phase. These are the type of fluctuations expected if the seeds for 
structure formation result from the superposition of quantum mechanical zero-point fluctuations 
of the scalar field that drove inflation of the early Universe, in the simplest inflation models, as 
discussed in Sec. [3] above. In the simplest inflation models the fluctuations are adiabatic. Fur- 
thermore, observational d ata are consistent with only adiabatic perturbations, so in what follows 
we focus on this case [see iBean et al.l (|2006l ) for a recent discussion of constraints on isocurvature 
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models] . 



As discussed in Sec. [3] above and Sec. 17.21 below, current large-scale observational results are 
reasonably w ell fit by an n = 1 scale-invariant pr imor dial spectrum of perturbations, the kind 
considered by iHarrisonl (|1970l ) , iPeebles <fe Yul (jl970l ) , and IZel'dovichl (119721 ) , and predicted in some 
of the simpler inflation models. The effect of causal physics on the later growth of structure, as 
discussed above, may then be represented by a "transfer function" that describes the relative growth 
of fluctuations on different wavelength scales. Observations of the anisotropy of the CMB and the 
clustering of galaxies and clusters at the present epoch probe the shape of the transfer function 
(as well as the primordial spectrum of perturbations) and thereby constrain structure formation 
models. Such observations are discussed below in Sees. 15.21 and [6l 



5.1.9. Gravitational Waves and Magnetic Fields 



As noted in Sec. [3] above, more complicated models of inflation can generate gravitational wave 
or magnetic field fluctuations that break adiabaticity. A primordial magnetic field might provide 
a way of explaining the origin of the uniform part of contemporary galactic magnetic fields; there 
are enough charged particles left over after recombination to ensure that pri mordial m a gnetic field 
lines will be pul led in, and the field amplified, by a collapsing gas cloud. iMaggiord (|2000l ) and 
Buonannd ( 12004) revie w primordia l grav ity waves, and cosmological magnetic fields are reviewed 
by IWidrowl (|2002l ) and iGiovanninil (|2004l ) . In the next subsection we consider the effects of such 
fields on the CMB. 



5.2. CMB Anisotropies 



As a result of the gravitational growth of inhomogeneities in the matter distribution, when the 
photons decouple from the baryons at last scattering at a redshift z ~ 10 3 (see Sec. 15.11 above) the 
photon temperature distribution is spatially anisotropic. In addition, in the presence of a CMB 
temperature quadrupole anisotropy, Thomson-Compton scattering of CMB photons off electrons 
prior to decoupling generates a linear polarization anisotropy of the CMB. After decoupling the 
CMB photons propagate almost freely, influenced only by gravitational perturbations and late-time 
reionization. Measurements of the temperature anisotropy and polarization anisotropy provide 
important constraints on many parameters of models of structure formation. This area of research 
has seen spectacular growth in the last decade or so, following the COB E discovery of the CMB 
temperature an i sotrop y. It has been th e subi ect of recent re views: see IWhite fe Cohnl (12 002') . 



(2005), and 



Ilu ,V Diu lelsoJ ^OoA IPeebles fc Eatral ibood . Sec. IV.B.ll), ISubramanianl (|2005h 
^Challinori (|2005h . Here we focus only on a few recent developments. 



Giovannini 



The three-year WMAP observations of CMB temperature anisotropies (jHinshaw et al.l 120071 ) 
are state-of-the-art data. On all but the very largest angular scales, the WMAP data are consistent 
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with the assumption that the CMB tempe rature anisotropy is well-described b y a spatial Gaus 



sian random p rocess (jKomatsu et al.ll2003l ). consistent with earlier indications ([Park et al.l 12001 



Wu et al.l l200ll ) . The few largest-scale angular modes exhibit a lack of power com pared to what 



is exp ected in a spatially-flat CDM model dominated by a cosmological constant ([Bennett et al 



2003al ). resulting in some debate about the assumptio ns of large-scale G aussianity and spatial 
isotropy. This feature was also seen in the COBE data (IGorski et al.lll998l ). The estimated large- 
angular-scale CMB temperature anisotropy power depends on the model used to remove foreground 
Galactic emission cont amination. Much work has been devoted to understanding foregro und emis- 
sion on all scales (e.g.. IMukherjee et alj|2003al : iBennett et alj|2003bl : iTegmark et al.ll2003l ). and the 
current consensus is that foregrounds are not the cause of the large- angular-scale WMAP effects. 

The CMB temperature anisotropy is conventionally expressed as an expansion in spherical 
harmonic multipoles on the sky, and for a Gaussian random process the multipole (or angular) 
power spectrum completely characterizes the CMB temperature anisotropy. The observed CMB 
anisotropy is reasonably well fit by assuming only adiabatic fluctuations with a scale-invariant 
power spectrum. These observational results are consistent with the predictions of the simplest 
inflation models, where quantum-mechanical fluctuations in a weakly-coupled scalar field are the 
adiabatic, Gaussian seeds for the observed CMB anisotropy and large-scale structure. 

Smaller-scale inhomogeneities in the coupled baryon-radiation fluid oscillate (see Sec. 15.11 
above), and at decoupling some of these modes will be at a maximum or at a minimum, giving rise 
to acoustic peaks and valleys in the CMB anisotropy angular spectrum. The relevant length scale is 
the acoustic Hubble length at the epoch of recombination; this may be predicted by linear physics 
and so provides a standard ruler on the sky. Through the angular diameter distance relation, the 
multipole numbers I of oscillatory features in the temperature anisotropy spectrum C( reflect space 
curvature (Ok) and the expansion history (which depends on Om and Oa) of the Universe. The 
angular scales of the peaks are sensitive to the value of the matter density parameter in an open 
Universe, but not in a spatially-flat (Ok = 0) Universe dominated by a cosmological constant, 
where the first peak is at a mult ipole index t ~ 220. This p rovid es a useful way to measure th e cur- 
vature of spatial hyper surf aces. ISugiyama k, Goudal (|1992l ) and iKamionkowski et al.1 (Il994albj) ar e 
early discussions of the C MB temperature anisotro py in an open Univers e , and iBrax et al.l (|2000l ). 
Baccigalupi et al. ( 2002 ). Caldwell Doran ( 2004 ). and Mukherjee et al. (j2003bl ) consider the case 
of scalar field dark energy in a spatially-flat Universe. CMB temperature an isotropy data on the 



position of the fir s t peak is consist ent with flat spatial hypersurfaces (e.g., iPodariu et al.1 12001 



Durrer et al.l 120031 ; iPage et al.l 120031). Model-based CMB data analysis is used to constrain mor e 
cosmological parameters (e.g., iLewis Bridld 120021 ; IMukherjee et al.l 12003a : ISpergel et al.l 120071 ) . 
For example, the relative amplitudes of peaks in this spectrum are sensitive to the mass densities 
of the different possible constituents of matter (e.g., CDM, baryons, and neutrinos, Ocdm> Ob, and 

The CMB pol arization anisotrop y was first detected from the ground by the DASI experiment 
at the South Pole (jKovac et al.ll2002l ). The three-year WMAP observations are the current state of 
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the art (jPage et al.ll2007l ). For a recent review of polarization measurements see lBalbi et al.l ((2006). 
The polarization anisotropy peaks at a larger angular scale than the temperature anisotropy, indi- 
cating that there are inhomogeneities on scales larger than the acoustic Hubble length at recom- 
bination, consistent with what is expected in the inflation scenario. The polarization anisotropy 
signal is interpreted as the signature of reionization of the Universe. The ability of WMAP to 
measure polarization anisotropies allows this experiment to probe the early epochs of non-linear 
structure formation, through sensitivity to the reionization optical depth r. 

Primordial gravitational waves or a primordial magnetic field can also generate CMB anisotropies 
Of particular current interest are their contributions to various CMB polarization anisotropies. (Be- 
cause polarization is caused by quadrupole fluctuations, these anisotropies constrain properties of 
the primordial fluctuations, such as the ratio of tensor to scalar fluctuations, r.) The effects of 
gravity waves on the CMB are disc ussed in the more recent standard cosmo logy and a s tropa rticle 
textbooks and bv iGiovanninil (120051 ). The magnetic field case is reviewed by Giovannini (2006) and 
Subramanian ( 20061): recent topics of interest m ay be traced from Lewis ( 2004 ). Kahniashvili &; Ratra 



(2005, 



20oJ, mid 



Brown fc Crittenden! (j2005l ). 



We continue discussion of the CMB anisotropies and cosmological parameters in Sec. 17.21 



5.3. Galaxy Formation and the End of the Dark Age 



The emission of the first light in the Universe, seen today as the CMB , is followed by a 
"dark age" before the first stars and quasars form. iBromm &: Larson! (12004! ) review formation 
of the first stars. Eventually, high energy p hotons from stars and quasars reionize interg a lactic 



2006 



Choudhurv k, Ferrara 



2006a 



Locb 



gas th roughout the Universe (for reviews see I Fan et al. 
2006a. b). Observations of polarization of microwave background photons by WMAP ([Page et al 
20071 ) suggest that reionization occurs at re dshift z ~ 11. However, strong absorption of Lyman-a 



redshift z ~ 6 (Becker et al 



2001 



photons by inte rgalactic neutral hydrogen (IGunn fc Peterson! Il965l ). seen in spectra of quasars at 



Fan et al 



2002T ) indicates that r eionization was not complete 



until somewhat later. This is an area o f ongoing research (see, e.g., IChoudhury Ferrara! l2006bl : 
Cnedin fc Fanlhood : I Alvarez et allhoOfih . 



Current models for galaxy formation follow the picture (jHovld Il953l : I Silk! Il977l ; iBinnev! Il977l ; 
Rees fc Ostrikell977i : Iwhite & Ree3l978h in which dark matter halos form by collisionless collapse, 
after which baryons fall into these potential wells, are heated to virial temperature, and then cool 
and condense at the centers of the halos to form galaxies as we know them. In short, baryons fall 
into the gravitational potentials of "halos" of dark matter, at the same time that those halos grow 
in size, hierarchically aggregating small clumps into larger ones. The baryons cool by emitting 
radiation and shed angular momentum, leading to concentrations of star formation and accretion 
onto supermassive black holes within the dark matter halos. 



In addition to the perturbative approach to structure formation discussed in Sec. 15.11 Lemaitre 
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also pioneered a "nonperturbative" approach ba sed on a spherically symmetric solution of the Ein 
stein equations. This spherical accretion model 



Gunn fc Gott [1972) describes the salient features 



of the growth of mass concentrations. See lGottJ (|1977l ). Peebles! ()1993l . Sec. 22), and lSahni Coles 
()1995l ) for reviews of such models. 



A phenomenological prescription for the statistics of non-linear collapse of str ucture, i.e., the 
form a tion of gravitationaily b ound objects, is given by the Press-Schechter formulae (jPress Schechter 
1974 ; ISheth &: Tormenl Il999l ) . Attempts to firm up the theoretical basis of such formulae form 
the "excursion set" forma lism which treats t h e formation of a gravitationaily bound halo as the 
result of a rand o m wa lk (|Mo fc White! Il996l ; IWhitel Il996l ; ISheth et al.1 l200ll ) . For a review see 
Coorav Shethl ((2002J) . These methods provide probability distributions for the number of bound 
objects as a function of mass t hreshold, and can be generalized to develop a complementary descrip- 
tion of the evolution of voids (ISheth van de Weygaertll2004l ). A more rigor o us approach assumes 



structure forms at high peaks in the smoothed density field (IKaiserl Il 984: 



Bar deen et al 
A vila- R eese! (hoOfih and 



1986; 



Baugh 



Sahni &; Colesl Il995l ) . Recent reviews of galaxy formation include 

(2006J). The next subsection, 15.41 describes numerical methods for studying structure formation. 



Apparent confirmation of the hierarchical picture of structure formation includes the striking 
images of galaxies appa rently in the process of assembly obtained by the HST in the celebrated 



"Hubble Deep Fields" ([Ferguson et al 



2000 



Beckwith et al.1 120061 ) . The detailed properties of 



galaxies and their evo l ution are o utside the s cope of this review. Recent reviews of the observat ional 
situation are iGawiserl (120061 ) and lEllisI (|2007l ) . Texts covering this topic include ISpinradl (|2005l ) and 



Longairi ( 200. 



While the current best model of structure formation, in which CDM dominates the matter 
density, works quite well on large scal es, current observa t ions indicate some poss ible problems with 
the CDM mode l on smaller scales; see iTasitsiorml (j2003l ) . [Peebles fc Ratral (hood . Sec. IV.A.2), and 
Primackl (|2005l ) for reviews. Sim ulations of structure formation indicate t hat CDM model halos 
may have cores that are cuspier (INavarro et al. 199 71: ISwaters et al.l l2000l ) and central densities 



that are higher (IMoore et al 



1999a 



Firmani et al. 



2001 



) than are observed in galaxies. Another 
concern i s that CDM models predict a larger tha n observed number of low- mass satellites of massive 
galaxies (jMoore et al.lll999bl ; iKlvpin et al.lll999l ). These issues have led to consideration of models 
with reduced small-scale power. However, it seems difficult to reconcile suppression of small-scale 
power with the observed small-scale clustering in the neutral hydrogen at redshifts near 3. 

The relationship between the distributions of galaxies (light) and matter is commonly referred 
to as "biasing." The currently- favored dark energy dominated CDM model does not require signifi- 
cant bias between galaxies and matter; in the best-fit model t he ratio of galaxy to matter clustering 
is close to unity for ordinary galaxies (jTegmark et al.l l2004al ) . 
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5.4. Simulations of Structure Formation 



Cosmological simulations usin g increasingly sophis ticated numerical methods provide a testbed 
for models of structure formation. iBertschingert (|1998l ) reviews methods and results. 



Computer simulations of structure formation in the Universe began with purely gravitational 
codes that directly compute the forces between a finite number of particles ("Particle-Particle" 
or PP codes) that sa mple the matter distribution. Early results used direct A-body calculations 
(jAarseth et al.l [l979j) . Binning the particles on a grid and computing the forces using the Fast 
Fourier Transform (the "Particle-Mesh" or PM method) is computationally more efficient, allowing 
simulation of larger volumes of space, but has force resolution of the order of the grid spacing. A 
compromise is the P 3 M method, which uses PM for large scale forces supplemented by direct PP 
calcul ations on small scales, as used for the important suite o f CDM simulations by 
(j 19851 ) . For details on these methods see iHocknev h Eastwood! (|1988l ) . 



Davis et al 



The force resolution of PM codes and the force resolution and speed of P 3 M codes may be 



1991 



does this dynamically t o increase force resolution in the PM gravity solver (IKravtsov et al. 
Norman Sz Brvanlll999l ). 



incre a sed by employing multiple gri d levels (|Villumsenlll989l ; lCouchmanlll991l;lBertschinger fc Gelb 
Gnedin Sz BertschingeJ 1996 ). Adaptive mesh refinement (AMR; Berger fc CollelaT 1989 ) 



1997 



Another approach to achieving both speed and go od force resolution in gravitational A^-body 
simulation is use of the hierarchical tree algorithm ( Barnes &: Hut! 19861). La rge cosmological 



simulations have used a parallelized version of this method (IZurek et al.l 1199^). Sig nificant in- 



crease in speed was found with the Tree Particle-Mesh algorithm (jBode et al.l l2000l ). GOTPM 
(jDubinski et al.l 12004 ). a parallelized hybrid PM+tree scheme , has been used for simulations in- 
volving up to 8.6 x 10 9 particles. PMFAST (|Merz et alJ 120051 ) is a recent parallelized multi-level 
PM code. 

Incorporation of hydrodynamics and radiative transfer in cosmological simulations has made it 
possible to study not only the gravitational formation of dark matter halos, but also the properties 
of baryonic matter, and thus the formation of galaxies associated with thos e halos. Methods for 
solving the fluid equations include smooth-particle hydrodynamics [SPH; see iMonaghanl (|l992l ) for 
a review], which is an inherently L agrangian app roac h, and Eulerian grid meth ods. Cosmological 
SPH simulations were pioneered by lEvrardl (|1988l ) and lHernquist h Katzl (|1989i ). To date, the cos- 
mological si mulation with the lar gest number of particles (10 10 ) employs SPH and a tree algorithm 
(GADGET; Sprineel et al.ll200ll). Grid-based codes used for cosmological simulation include that 



described bylCen] (119921 ) and lPvu et al.1 (119931 ) 



To date, no code has sufficient dynamic range to compute both the large scale cosmological 
evolution on scales of many hundreds of megaparsecs and the formation of stars from baryons, but 
physical heuristics have be en success fully incorporated into some codes to model the conversion of 
baryons to stars (see, e.g.. ICenlll992l ). 
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The Millenium Run simulation (jSpringel et al.ll2005l ) represents the current state-of-the art in 
following the evolution of both the dark matter and baryonic components on scales from the box 
size, 500/i -1 Mpc, down to the resolution limit of roughly 5/i _1 kpc. See this article and references 
therein for discussion of the many pieces of uncertain physics necessary for producing the observed 
baryonic structures. 

Another approach to modeling the properties of the galaxies associated with dark matter ha- 
los is to use the hi s tory of halo mergers together with semi-analytic modeling of galaxy p roperties 
(|Lacev et al.lll993l : Ecauffmann et aDll993l : Icole et al.lll994l : ISomerville k PrimacM Il999l ). When 
normalized to the observed luminosity function of galaxies and Tully-Fisher relation for spiral 
galaxies, these semi-analytic models (SAMs) reproduce many of the observed features of the galaxy 
distribution. A common approach is to use SAMs to "paint on" the properties of galaxies that 
would resid e in th e dark matter halos found in purely gravitational simulations. See 



Avila-Reese 



(J2006J) and iBaughl (120061) for recent reviews. Related t o the SAMs approach are halo occupation 



models (jBerlind Weinberg 



2002 



Kravtsov et alJl2004l ) that parameterize the statistical relation- 



ship between the masses of dark matter halos and the number of galaxies resident in each halo. 



6. Mapping the Universe 

The observed features of the large-scale distribution of matter include clusters, superclusters, 
filaments, and voids. By mapping the distribution of galaxies in the Universe, both in two dimen- 
sions as projected on the sky and in three dimensions using spectroscopic redshifts, astronomers 
seek to quantify these inhomogeneities in order to test models for the formation of structure in 
the Universe. Not only the spatial distribution of galaxies, but also the distribution of clusters of 
galaxies, quasars, and absorption line systems provide constraints on these models. Peculiar veloci- 
ties of galaxies, which reflect inhomogeneities in the mass distribution, provide further constraints. 
Here we briefly review important milestones and surveys relevant for testing cosmological models. 



6.1. Galaxy Photometric Surveys 



Studies of the global spacetime of the Universe assume the "cosmological principle" ((Milne 



19331 ) which is the supposition that the Universe is statistically homogeneous when viewed on 
sufficiently large scales. The angular distribution of radio galaxies provides a good test of this 
approach to homogeneity, because radio-bright galaxies and quasars may be seen in flux-limited 
samples to nearly a Hub ble distance, c/Hq. Indeed , the ~ 31,000 brightest radio galaxies observed 
at a wavelength of 6 cm ((Gregory Condon 119911) are distribu ted nearly isotropically, and similar 
results are found in the FIRST radio sur yev (IBecker et alJll995T ). [For a review of other evidence for 
large-scale spatial isotropy see Sec. 3 of lPeeblesI Jl993l ).] In contrast, the Universe is clearly inho- 
mogeneous on the more modest scales probed by optically-selected samples of bright galaxies, For 
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example, si gnificant clustering is observed among the roughly 30,000 galaxies in the lZwickv et al 
(|l961-1968h catalog. 



Maps of the distribution of nebulae revealed anisotrop y in the sky b efore a stronomers came 
to agr ee that many of these nebulae were distant galaxies (jCharlierl Il925l ) . The IShapley fc Ames 
(|1932l ) catalog of galaxies clearly showed the nearby Virgo cluster of galaxies. Surveys of selected 
areas on the sky using photographic plates to detect distant gala xies clearly revea l ed anisotropy of 
the ga laxy distribution and were used to quantify this anisotropy (|Mowbravlll938l ). Ide Vaucouleurs 
(|l953l ) recognized in this anisotropy the projected distribution of the local supercluster of galaxies. 



Rubinl ()1954l ) used two-point correlations of galaxy counts fro m Harvard College Observ atory 
plates to detect fluctuations on the scale of clusters of galaxies. The lShane Wirtanenl (|1954l ) Lick 
Survey of galaxies used counts of galaxies found on large-format photographic plates taken at Lick 
Observatory to make the first large-scale map of the angular distribution of galaxies suitable for 
statistical analysis. Early analysis of th ese data inclu d ed methods such as cou nts-in-cell analyses 
and the two-point correlatio n function (ILimberl 19541 ; iTotsuii fc Kiharal Il969l ). The sky map of 
the Lick counts produced by ISeldner et al.l (|l9771 ) visually demonstrated the rich structure in the 
galaxy distributio n. Peebles and collabor ators used these data for much o f their extensive work on 
galaxy clustering ( Groth fc Peebles 1977 ): for a review see Peebles (jl98ol . Sec. III). 



The first Palomar Observato ry Sky Survey (POSS) yie lded two important catalogs: the lAbell 



(|l958l ) catalog of clusters and the IZwicky et al.l (11961-19681 ) catalog of clusters and galaxies identi- 
fied by eye from the photographic plates. lAbelll (| 19611 ) found evide nce for angular "super clus tering" 
(clustering of galaxy clusters) that was confirmed statistically by lHauser Peebles! (|1973l ) . Pho- 
tographic plates taken at the UK Schmidt telescope were digitized using the Automatic Plate 
Measuring (APM) machine to produce the APM catalog of roughly two million galaxies. Calibra- 
tion with CCD photometry made the APM cat alog invaluable for accurate study of the angular 
correlation function of galaxies on large scales (jMaddox et al.lll990l ). Perhaps the last large-area 
galaxy photometric su rvey to employ photographic plates is the Digitized Palomar Observatory 
Sky Survey (DPOSS) (|Gal et alJbooj i. 



The largest imaging survey that employ s a camera with arra ys of charge-coupled devices 
(CCDs) is the Sloan Digital Sky Survey (SDSS: IStoughton et al.ll2002i ). The imaging portion of this 
survey includes five-color digital photometry of 8000 deg 2 of sky, with 215 million detected objects. 
Imaging for the SPSS is o btained using a special-purpose 2.5 m telescope with a three-degree field 
of view (IGunn et al.ll2006l ). 



Important complements to optical surveys include large-area catalogs of galaxies selected in 
the infrared and ultraviolet. Nearly all-sky source catalogs were produced from in frared photometry 
obtained with the Infrared Astronomical Sat ellite (IRAS; Beichman et al.lll988l ) and the ground- 
based Two Micron All Sky Survey (2MASS: barrett et al.l I2OO0I ) . The ongoing Galaxy Evolution 
Explorer satellite (GALEX; iMartin et al.ll2005l ) is obtaining ultraviolet imaging over the whole sky. 
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6.2. Galaxy Spectroscopic Surveys 



Systematic surveys of galaxies using spectroscopic redshifts to infer their distances b e gan with 
obser vations of galaxies selected from the Shapley-Ames catalog (jHumason et al.l Il956l ; ISandage 
19781 ). Important early mapping efforts includ e identification of su perclu sters and voids in the 



distrib ution of galaxies and Abell clusters by iJoeveer et al.l (119781 ) . the iGregory &: Thompson 



(197 8]) study o f the Coma/Abelll367 supercluster and its environs that identified voids, and the 
Kirshner et al.l (|198ll ) study of the correlat ion function of galaxies and discovery of the giant void in 
Bootes. Early targeted surveys include the lGiovanelli fc Haynea (| 19851 ) survey of the Perseus-Pisces 
supercluster. 

Redshift surveys of large are as of the sky began with the first Center for Astrophysics redshift 
survey (CfAli lHuchra et alJll983l ). which includes redshifts for 2401 galaxies brigh ter than apparent 
magnitude = 14.5 over a wide area toward the North Galactic Pole. CfA2 (jFalco et al.l 11999) 
followed over roughly the same area, extending to apparent magnitude = 15.5. At this depth. 



the r ic h pattern of voids, clusters , and superclusters were strikingly obvious (Ide Lapparent et al 



19861 ). iGiovanelli Haynesl (Il99ll ) review the status of galaxy redshift surveys ca. 1991. 



Both CfA redshift surveys used the Zwi cky catalog of galaxie s to select targets for spectroscopy. 
The Southern Sky Redshift Survey (SSRS; lda Costa et alJll99a ) covers a large area of the south- 
ern hemisphere (contiguous with CfA2 in the northern galactic cap) to similar depth, using the 
ESO/Uppsala Survey to select galaxy targets and a spectrograph similar to that employed for the 
CfA surveys. The Optical Redshift Survey (ORS) supplemented existing redshift ca talogs with 
1300 new spectroscopic redshifts to create a nearly all-sky survey ([Santiago et alJll995l ). 



Deep "pencil-beam" surve ys of narrow patches on the sky revealed apparently-periodic struc- 
ture in the galaxy distribution (jBroadhurst et alJll99d ). 



The Las Campanas Redshift Survey (LCRS: IShectman et alJll99a ). the first large-area survey 
to use fiber optics, covered over 700 deg 2 near the South Galactic Pole. This survey was important 
because it showed that structures such as voids and superclusters found in shallower surveys are 
ubiquito us but the structu res seen by LCRS were no larger th an those found earlier . The Century 
Survey (jGeller et al.l Il997h and the ESO Deep Slice survey (jVettolani et al.l Il998f ) were likewise 
useful for statistically confirming this emerging picture of large-scale structure. 

Sparse surveys of galaxies to efficien tly study statistical properties of the galaxy distribution 
include the Stromlo-APM redshift survey (jLoveday et alJll99o ) based on 1/20 samp ling of the APM 
galaxy catalog and the Durham/UK Schmidt redshift survey (jRatcliffe et al.lll998l ). 



While optically-selected surveys are relatively blind to structure behind the Milky Way, redshift 
surveys based on objects detected in the infrared provide nearly all-sky coverage. A sequen ce of 
sur veys of objects dete cted by IRAS were carried out, flux -limited to 2 Jy (jStrauss et alJll992n. 1.2 



Jv (jFisher et al.lll995l ). and 0.6 Jy (ISaunders et alJl2000l ). The 6dF Galaxy Survey (IJones et al 



20041 ) will measure redshifts of 150, 000 galaxies photometrically identified by 2MASS (IJarrett et al 
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2000). 



The 2-degree Field Galaxy Redshift Survey (2dFGRS) of 250, 000 galaxies (IColless et al 



2001 



was selected from the APM galaxy catalog and observed using the two-degree field multi-fiber 
spectrograph at the Anglo- Australian 4 m telescope. The survey is complete to apparent magnitude 
mj = 19.45 and covers about 1500 deg 2 . 



The spectroscopic component of the SDSS (jStoughton et al.ll2002l ) includes medium-resolution 
spectroscopy of 675, 000 galaxies and 96, 000 quasars identified from SDSS photometry. These 
spectra are obtained with dual fiber-optic CCD spectrographs on the same 2.5 m telescope. The 
main galaxy redshift survey is complete to m r = 17.77 and is complemented by a deeper survey of 
luminous red galaxies. The ongoing extension of this survey (SDSS-II) will expand the spectroscopic 
samples to more than 900, 000 galaxies and 128, 000 quasars. 

Spectroscopic surveys t hat trace structu re in the gal axy distribution at much larger redshift 
include the DEEP2 survey (jCoil et al.1 12004 ) an d others (ISteidel et all 12004 ) employing the Keck 
Observatory, and the VIMOS VLT Deep survey (|Le Fevre et al.l l2005). 



6.3. Cluster Surveys 



Mapping of t he Universe using galaxy clusters as tracers began with study of the Abell cat- 
alog (lAbelll Il958l; lAbell et al.lll989l ). Studies of the angular clustering of Abell clusters includes 
Hauser fc Peebles! (119731). Several red s hift s urvey s of Abe l l clus ters have been conducted, includ- 
ing those described bv lPostman et al.1 (|1992l ) and iKatgert] (|1996l ). Important cluster samples have 
also been identified from digitized photog raphic plates from the UK Schmidt tele scope, followed 
up by redshift surveys of cluster galaxies (jLumsden et al.l Il992l ; I D alt on et al.l I1992T ) . More distant 



samples of clusters have b een identified using deep CCD photometry (see, e.g., iPostman et al 



1996 



Gladders Yeell2005l ). In X-ray bandpasses, clus ter samples useful for studying large-scal e 



structure have been identified using data from ROSAT (IRomer et al.lll994l : iBohringer et al.ll2004h . 
The SDSS is produ cing large catalogs of galaxy clusters using a variety of selection methods 



( Bah call et al 



20031 ). Use of the Sunyaev-Zel'dovich effect (the microwave decrement caused by 



Thomson-Compton scattering of the CMB p hotons by the intraclu ster gas) holds great promise to 
identify new dee p samples of galaxy clusters (ICarlstrom et alj|2002f>. General reviews of clusters of 
galaxies include iRosati et al.l (j2002j) , IVoitl (|2005l ) , and iBorganil (|2006l ) . 



6.4. Quasar surveys 



The advent of multi-object wide-field spectrographs has made possible collection of very large 
samples of spectrq scopically-confir med quasars, as obser ved by the 2dF QSO Redshift Survey 
(|Croom et al.ll2004h and the SDSS (jSchneider et allboosl ). For a ca. 1990 review of the field see 
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Hartwick k Shadd (|1990l ). While quasars themselves are too sparsely distributed to provide good 
maps of the l arge-scale distribution of matter, their clustering in red shift space has been measured 
(|Osmerlll98ll ) and generally found to be similar to that of galaxies (jOutram et al.l 120031) . Similar 



resul ts obtain from clustering analyses of active galactic nuclei in the nearby universe (IWake et al 



2004 ). although this clustering depends in detail on the type of AGN (jConstantin k Vogelevll2006l ) 



The distribution of absorption lines from gas, particularly fr om the Lyman-q "forest" of neutral 
hydrogen clouds along the line of sight toward bright quasars (Lynds 1971 ; Rauch 199 prov ides 
an important statistical probe of the distribution of matter (see, e.g., iMcDonald et al.ll2005l ) on 
small scales and at large redshift. 



6.5. Peculiar Velocity Surveys 



When measured over sufficiently large scales, the peculiar motions of gal axies or cluste r s sim 
ply depend on the p otential field generated by the mass distribution (see iPeebled Il980i . 



1993 



m 



Davis k Peebles! 1 19831). Te c hniques for measurin g distances to other galaxi es are critica l ly rey iewed 



Eowan-Eobinsonl 111983 ). Uacobv et al.1 (|1992 ). IStrauss k Willicki (|l995l ). and lWebbl ((1999). To 



gether with the galaxy or cluster redshifts, these measurements yield maps of the line-of-sight 
component of the pec uliar velocity. From such da t a it is possib le to reconstruct a map of the mat- 



ter density fi eld (e.g. , Bcrtsching er k Deke] 



in time (e.g., iPeebles 



1990 



1989 



Dekell Il994l ) or to trace the galaxy orbits back 



Goldberg k Sper gcl 2000). Analyses of cor relations of the den sity and 



velocity fields also yield constraints on the cosmic matter density (e.g.. IWillick et al.lll997l ). 



Dressier et al 



Rubin et a .1 (119761 ) were the first to find evidence for bulk flows from galaxy peculiar velocities. 

(fl98?T ) f ound evidence for a bulk flow toward a large mass concentration, dubbed 



the "Great Attractor." lLauer k Postman! ()1994l ) found surprising evidence for motion of the Local 
Group on a larger scale. However, analysis of subsequent peculiar velocity surveys indicates that 
the inferre d bulk flow results, i ncluding those of Lauer and Postman, are consi stent within the un 



certai nties (jHudson et al.ll2000n. The status of thi s field ca. 1999 is su rvey ed bvlCourteau k Wi llick 



d200fl), recent results include iHudson et al.1 (|200J), and[Dekej 4l99J) and lStrauss fc Willicki (1199m 
review this topic. Comparison of peculiar velocity surveys with the peculiar velocity of our Galaxy 
implied by the CMB dipole indicates that a significant component of our mo tion must arise from 
density inhomogeneities that lie at rather large distance, beyond 60/i _1 Mpc ( Hudson et al.1 12004 ). 



6.6. Statistics of Large-Scale Structure 

The clustering pattern of extragalactic objects reflects both the initial conditions for structure 
formation and the complex astrophysics of formation and evolution of these objects. In the standard 
picture described above, linear perturbation theory accurately describes the early evolution of 
structure, thus measurement of clustering on very large scales, where the clustering remains weak, 
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closely reflects the initial conditions. On these scales the density field is very nearly Gaussian 
random phase, therefore the two-point correlation function of the galaxy number density field (also 
called the autocorrelation or covariance function) or its Fourier transform, the power spectrum, 
provides a complete statistical description. (Temperature anisotropies of the CMB discussed in 
Sec. 15.21 arise from density fluctuations at redshift z ~ 10 3 that evolve in the fully linear regime.) 
On the scales of galaxies and clusters of galaxies, gravitational evolution is highly non-linear and the 
apparent clustering depends strongly on the detailed relationship between mass and light in galaxies. 
In between the linear and non-linear regimes lies the "quasi-linear" regime in which clustering 
growth proceeds most rapidly. A wide range of statistical methods have been developed to quantify 
this complex behavior. Statistical properties of the gala xy distr i bution and details of est i matin g 
mos t of the relevant statistic s are described in depth by Peebles! (|1980l ). iMartmez Saarl ((2002), 



and Bernardeau et al 



(2002 V Meth ods of using galaxy redshift surveys to constrain cosmology 



are reviewed bv lLahav Sutol (|2004l ) and iPercivall (|2006l ). Constraints on cosmological parameters 
from such measurements are discussed below in Sec. 17.21 

The simplest set of statistical measures are the n-point correlation functions, which describe 
the joint probability in excess of random of finding n galax ies at specified relative separation . Earl y 
app lications of correla t ion fu nctions to galaxy data include lLimberl (|1954l ). lTotsuii Kiharal (|1969l ). 
and iGroth &; Peebles! (|1977l ). The n-point functions may be estimated by directly examining the 
positions of n-tuples of galaxies or by using moments of galaxy counts in cells of varying size. Test s 
of scaling relations among the n-point functions are discussed in detail bv lBernardeau et al.l ((2002). 



Power spectrum analys es of large galaxy redshift surveys (jVogeley et al 



1993 



1992 



Fisher et al 



Tegmark et al.l l2004bl ) yield useful constraints on cosmological models. Closely related to 



power spectrum analys e s are estimates o f cosmological p a ramet ers using orthogonal functions 



(jVogeley fc Szalav 



1996 



Verde et al 



Pope et al.l 12004) . [T egmark et al.l (|1998l ) discuss the merits of different 



(|2002l ) describe a measurement of the galaxy 



methods of power spectrum estimation, 
bispectrum. 

A number of statistics have been developed to quantify the geometry and topology of large- 
scale structure. The topological g enus of isodensity contours characterizes the connectivity of large- 



scale structure (Gottetal 



19871 ). Measurements of the genu s are consistent w ith random phase 



initial conditions (as predicted by inflation) on large scales (jGott et al.l [l989j) , with departures 
from Gaussi anity on smaller scales where nonlinea r gravitational evolution and biasing of galaxies 



Gaussianity of the CMB anisotropy ((Park et al 



are evident (jVogeley et al.l Il994 ; iGott et al.l 12006 \. Si milar techniques are used to check on the 



2001 



Wu et al 



2001 



well as identify foreground emission signals in CMB anisotropy data (IPark et al 



Komatsu et al 
2002). 



2003) 



as 



T he void prob ability function, which characterizes the frequency of completely empty regions of 
space dwhitelll979h. ha s been estimated from galaxy redshift surveys (jMaurogordato fc Lachieze-Revl 



1987 



algorithms (|E1-Ad et al 



Hovle &; Vogelevl 1 20041) . Catalogs of voids have been constructed with objective void finding 



1996 



Hovle fc Vogelevl 12002 ). 
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Early investigations of the pattern of galaxy clustering dating back to ICharlierl (119251) sug- 
gested a clustering hierarchy. The fractal model of clustering introduced by iMandelbrotl (|1982l . 
and references therein) further motivated investigation of the possibility of scale-invariant cluster- 
ing of galaxies. R e sults of suc h analys e s of g alaxy survey data were con troversial [compare, e.g., 
Svlos Labini et al. I (jl998h with iHattonl (|1999l ) and iMartmez et al.l (|200ll ) and references therein] . 



While fractal behavior is seen on small scales, there is fairly strong evidence for an approach 
to homogeneity in galaxy redshift and photometric surveys on very large scales. Thus, a simple 
scale-invariant fractal description seems to be ruled out. A multi- fractal description of clu stering 
continues to provide a useful complement to other statistical descriptors ([ Jones et al.ll2005l ). Con- 
sideration of modified forms of the fractal picture are of interest for providing slight non-Gaussianity 
on very large scales that might be needed to explain the very largest structures in the Universe. 

Anisotropy of galaxy clustering in redshift space results from bulk flows on large scales that 
amplify clustering along the line of sight to the observer and from motions of gal axies in viria l- 
ized systems su ch as clusters that elongate those structures along th e line of sight (|Kaiserlll987l ). 
Hamilton! (|1998l ) provides an extensive review and lTinker et al.l ([2006) describe recent methods for 
estimating cosmological parameters from redshift-space distortions of the correlation function or 
power spectrum. 

The dependence of clustering statistics on properties of galaxies provides important clues to 
their history of formation and reflects the complex relationship between the distributions of mass 
and l uminous matter. Th e amplitude of gala xy clustering is seen to va ry with galaxy morphology 



(e.g.JDavis Gellerlll976l : lGuzzo et al.lll997l ) and with luminosity (e.g. JHamirtonlll988l : IPark et al 



19941 ). In recent analyses of the SPSS and 2dFGRS, these and similar trends with color, surface 



brightness, and spectral type are seen (jNorberg et al 



2002 



Zehavi et al.1 120051 ). 



Spectroscopy obtained with 8-10 m class telescopes has rec ently made it pos s ible to accurately 



study structure in the ga laxy distribution at higher redshift (jCoil et al 



2005 



Le Fevre et al.l 1200.51 ). 



2004 



Adelberger et al 



7. Measuring Cosmological Parameters 

7.1. The Case for a Flat, Accelerating Universe 

As mentioned in Sec. HI observations of Type la supernovae (SNela) provide strong evidence 
that the expansion of the Universe is accelerating. Type la supernovae have the useful property 
that their peak intrinsic luminosities are correlated with how fast they dim, which allows them to be 
turned into standard candles. At redshifts approaching unity, observations indicate t hat they are 



dimmer (and so farther away) than would be predicted in an unaccelerating Universe (jRiess et al 



1998 



Perlmutter et al.lll999l ). In the context of general relativity this acceleration is attributed 



to dark energy that varies slowly with time and space, if at all. A mass-energy component that 
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maintains constant (or nearly constant) density has negative pressure. Because pressure contributes 
to the active gravitational mass density, negative pressure, if large enough, can overwhelm the 
attraction caused by the usual (including dark) matter mas s density and resu lt in accelerated 
expansion. For a careful review of the early s upernova tests seelLeibundgutl (120011). For discussion s 
of the cosmological implications of this test see lPeebles fc Ratral (|2003l ) and iPer ivolarop oulod (|2006l ). 
Current supernova data show that models with vanishing cosmological constant are more than four 
standard deviations away from the best fit. 

The supernova test assumes general relativity and probes the geometry of spacetime. The 
result is confirmed by a test using the CMB anisotropy that must in addition assume the CDM 
model for structure formation discussed in Sec. H] (see Sec. 15.31 for apparent problems with this 
model). As discussed in Sec. 15. 2[ CMB anisotropy data on the position of the first peak in the 
angular power spectrum are consistent with the curvature of spatial hypersurfaces being small. 
Many independent lines of evidence indicate that the mass density of nonrelativistic matter (CDM 
and baryons) — a number also based on the CDM structure formation model — is about 25 or 30 
% of the critical Einstein-de Sitter density (see Sees. EH and \7.2h . Because the contemporary mass 
density of radiation and other relativistic matter is small, a cosmological constant or dark energy 
must contribute 70 or 75 % of the current m ass budget of the Univ erse. For reviews o f the C MB 



data constraints see iPeebles fc Ratral (|2003l ) , ICopeland et al.l (|2006l ) , and ISpergel et al.l (|2007l ) . 



7.2. Observational Constraints on Cosmological Parameters 



The model suggested by the SNela and CMB data, spatially flat and with contemporary mass- 
energy budget split between a cosmological constant or dark energy (~ 70 %), dark matter (~ 25 
%), and baryonic matter (~ 5 %), is broadly consistent with the results of a large number of other 
cosmological tests. In this subsection we present a very brief discussion of some of these tests and 
the constraints they impose on the parameters of this " stand ard" cosmologi c al mo d el. Two nice 
reviews of the cosmological tests are Sec. 13 of IPeebles! (|1993l ) and ISandage I (|1995l ). iHoga (|1999i ) 
provides a conci s e sum mary of various geometrical measures used in these tests. Section IV of 
Peebles fc Ratral (|2003l ) reviews more recent developments and observational constraints. Here we 
summarize some of these as well as the si gnificant progress of the last four years. Numerical values 
for cosmological parameters are listed in lLahav &; Liddld (|2006l ) , although in some cases there is 
still significant ongoing debate. 



There have been many — around 500 — measurements of the Hubble constant Hq, (jHuchra 



20071 ). the current expansion rate. Since there is debate about the error estimates of some of 
these mea surements, a med ian statistics meta-analysis estimate of Hq is probably the most robust 
estimate (|Gott et alJboOlh . At two standard deviations this gives H = WOh km s" 1 Mpc" 1 
= 68 ± 7 km s" 1 Mpc -1 = (14 ± 1 Gyr) -1 (|Chen et al.ll2003l ). where the first equat ion defines h. It 



is sig nificant that this result agrees with the estimate from the HST Key Project (jFreedman et al 



2001 



), the HST estimate of Sandage and collaborators (ISandage et al.ll2006l ). and the WMAP three- 
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year data estimate (which assumes the CDM structure formation model) (jSpergel et al.l 120071 ) 



A measurement of th e redshift dependence of the Hubble para meter can be used to constrain 



cosmological parameters (Jimene z Sz Loeb 



2002; 



using preliminary data see ISamushia &: Ratral (120061 ) and ISen Scherrerl (|2007l ) . 



Simon et al 



2005) . For applic ations of this test 



Expansion time tests are reviewed in iPeebles fc Ratral ([20031 . Sec. IV. B. 3). A recent devel- 
opment is the WMAP CMB an isotropy data estim ate of the age of the Universe, to = 13.7 ± 0.3 
Gyr at two standard deviations (jSpergel et al.ll2007l ). which assumes the CDM structure formation 
model . This WMAP to estima t e is consistent with to estimated from g lobular cluster observa- 



tions (IKrauss Chabovei 



2003 



Gra tton et al.l 120031 : llmbriani et al.l 12004 ) and from white dwarf 



star measurements ([Hansen et al.l 120041 ) . The above values of Hq and to are consistent with a 
spatially-flat, dark energy dominated Universe. 



As discussed in Sec. [231 IPeebles & Ratral (|2003 ) . iField Sarkarl (j2006T ). and ISteigmanl (|200J), 
4 He and 7 Li abundance measurements favor a higher baryon density than the D abundance measure- 
ments and the WMAP CMB anisotropy data. (This difference is under active debate.) However, 
it is remarkable that high-redshift (z ~ 10 3 ) CMB data and low-redshift (z £ few) abundance 
measurements indicate a very similar baryon density. A summary range of the baryonic den- 
sity parameter from nucleosynthesis is Qb = (0.0205 ± 0.0035)/i~ 2 at two standard deviations 
([Field fc Sarkaribood ). 



As mentioned above, Type la supernovae ap parent magnitude s as a function of redshift may 
be used to constra in the cosm o logica l mo del. See IPeebles k, Ratral ([20031 . Sec. IV. B. 4) for a sum- 



mary of this test. iRiess et al.l (| 19981 ) and IPerlmutter et al 



a cosmological constant from this test, and iPodariu &: Ratral (|2000l ) and IWaga Friemanl ([2000 ) 



1999) provi ded initial constra i nts on 



gener alize the method t o con st rain scalar field dark en e rgy. Developments may be trac e d bac k 
from IWang fc Tegmark! (hoojUciocAiatti et alJ ifooOfih . Ustier et al.1 iboOfjhJjfass et al.1 JiooJ), 



Nesseris Sz Perivolaropoulosl ( 20051 ). Jassal et al. ( 20od ). and iBarger et al.l fcooj ). Proposed satel 
lite experiments are u nder active discussion and should result i n tigh t cons traints on dark energy 
and its evolution. See IPodariu et al.l (|200ll ). IPerlmutter et al.l (|2006l ) . and iRefregier et al.l ([20061 ) 
for developments in this area. 

The angular size of objects (e.g., quasars, compact radio sources, radio galaxies) as a function 
of redshift provides another cosmological test. These observations are not as numerous as the 
supernovae, so this test is much less constraining, but the results are consist ent with those f r om the 
SNI a apparent magnitud e t est. Developments may be traced back through IChen Ratral (|2003al ) 



and IPodariu et al.l ([20031 ). iDalv &: Diorgovskil (j2006l ) describe a way of combining the apparent 



magnitude and angular size data to more tightly constrain cosmological parameters. 

"Strong" gravitational lensing, by a foreground galaxy or cluster of galaxies, produces mul- 
tiple images of a background radio source. The sta tistic s of strong len sing may be used to con- 
strain the cosmological model. iFukugita et al.l ()1990l ) and I Turner! ([1990] ) note that for low nonrel- 
ativistic matter density the predicted lensing rate is significantly larger in a cosmological constant 
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dominate d spatially-flat m o del th an in an open model. The scalar field dark energy case is dis- 
cussed in iRatra &: Quilled (|1992l ) and lie s between t hese two limits. For reviews of the test see 
Peebles fc Ratral (120031. Sec. IV.B.6) a nd iKochanekl (|2006l ). Recent developments may be traced 
back from iFedeli &: Bartelmannl (120071) . Cosmqlogical con strain ts from the CLASS gravitational 
lens statistics data are discussed in lChae et al.l (|2002l . 12004 ) . and lAlcaniz et al.l (|2005l ). These con- 
straints are consistent with those derived from the supernova apparent magnitude data, but are 
not as restrictive. 

Galaxy motions respond to fluctuations in the gravitational potential, thus peculiar velocities 
of galaxies may be used to es timate the non relat ivistic matter density pa rameter Qm [as discussed in 
Sees. [Hand [63J above and in lPeeblesI (|l999l ) and Peebles k, Ratral (|2003l . Sec. IV.B.7)] by comparing 
the pattern of flows with maps of the galaxy distribution. Note that peculiar velocities are not 
sensitive to a homoge neously distributed ma ss-energy component. For a summary of recent results 
from the literature see lPike Hudson! (|2005l ). Measurements of the anisotropy of the redshift-space 
galaxy distribution that is pr oduced by peculiar velocities also yield estimates of the matter density 



( Kaiser 



1987 



Hamiltonlll998l ). Most methods mea sure this anisotropy in the galaxy autocorr e lation 



or p ower spectrum (s e e, e.g ., iTinker et al.ll2006l ). Recent analyses include lHawkins et al.l ([2003) 



and Ida Angela et al. I (|2005l ) from the 2dFGRS and 2QZ surveys. Also of interest are clustering 
analyses of the SDSS that explicitly take into account t his redshift-space anisotropy either by using 



the predicted distortions when constructing eigenmod es dPope et al 



that are sensitive to radial vs. angular fluctuations ([Tegmark et al 



2004) or by constructing modes 



2004b|). 



A median statistics analysis of density estimates from peculiar velocity measurements and a 
variety of other data indicates that the nonr elativistic matter dens ity parameter lies in the range 



0.2 ^ ~ 0.35 at two standard deviations ([Chen &: Ratral l2003bl ) . This is consist ent with esti- 



mate s based on other data, e.g., the WMAP CMB data result in a very similar range ([Spergel et al 
20071 ). 



"Weak" gravitational lensing (which mildly distorts the images of background objects), in 
combination with other data, should soon provide tight constraints on t he nonrela t ivistic mat- 
ter density paramet er. F or reviews of weak l e nsing see iRefregierl ( Eooah ISchneiderj (hood), and 
Munshi et all feood ). See ISchimd et all (l2007h . iHetterscheidt et all l|2007l ). iKitching et aD (|2007h 
for r ecent developmen t s. Weak gravitatio nal lensing also provides evidence for dark matter (see, 



e.g. 



Clowe et al 



2006 



Massev et al.ll2007h . 



Rich clusters of galaxies are thought to have originated from volumes large enough to have fairly 
sampled both the baryons and the dark matter. In conjunction with the nucleosynthesis estimate 
of the baryonic mass density parameter, the rich cluster estimate of the ratio of baryonic and 
nonrelativistic (including baryonic) matter — the cluster baryon fraction — provides an estimate 



of the nonrelativistic matter density parameter (jWhite et al.lll993l ; iFukugita et al.lll998l ) . Estimates 
of from this test are in the range listed above. A promisi ng method for measuri ng the cluster 



baryonic gas mass fraction uses the Sunyaev-Zel'dovich effect (jCarlstrom et al.ll2002l ). 
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An extension of this cluster test makes use of measurements of the rich cluster baryon mass 
fraction as a function of redshift. For relaxed rich clusters (not those in the process of collapsing) 

The cluster baryon fraction depends on 



( Sasaki 


1996; 


Pen 


1997) 



clusters at the right angular diameter distances to ensure that the cluster baryon mass fraction is 
independent of redshift. This test provides a fairly restri ctive constra i nt on Qm, consistent with th e 



range above; developm ents may be traced bac k through I Allen et al.l (120041 ) , IChen k Ratral (J200J) 



Kravtsov et al.l (|2005l ). and lChang et al.1 (120061 ). When combined with complementary cosmological 
data, especially the restrictive SNIa data, the cluster baryon mass fraction versus redshift data 
provide tight constraints on the cosm ological model, favoring a cosmological constant but not ye t 



ruling out slowly varying dark energy (jRapetti et al. 



2005 



dl2 



Alcaniz k Z hu 2005: 1 Wilson et al.ll2006l ') . 



The number density of rich clusters of galaxies as a function of cluster mass, both at the 
present epoch and as a function of redshift, constrains the a mplitude of mass fluct uations and the 
nonrelativistic matter density parameter (see Sec. IV. B. 9 of lPeebles k Ratral 120031 . and references 
therein) . Current cluster da t a favor a matter density paramet er in the range discussed above 



(Rosatietal 



2002 



Voitll2005l : lYounger et al. 



2005 



Borganj|2006). 



The rate at which large-scale structure forms could eventually provide another direct test 
of the cosmological model. The cosmological constant model is discussed in iPeeblesI ()1984j) and 
some of the more recent textbooks listed b elow. The sca l ar fie ld dark energy model is not a s 
trac table; developmen ts may be traced from iMainini et al.1 (|2003l ) , iMota k van de Bruckl (J2004J) , 
and iMaio et al.l (hood ) . 



Meas urements of CMB temp erature and polarization anisotropics [see Sec. 15.21 above and Sec. 
IV.B.ll of lPeebles k Ratral ([2003)] provide some of the strongest constraints on several cosmological 
model parameters. These constraints depend on the assumed structure formation model. Current 
constraints are usually base d on the CDM mode l (or some variant of it). As discuss ed in Sec . 15.2 



the t hree-year WMAP data (jHinshaw et al.ll2007l ) provide state-of-the-art constraints (jSpergel et al 
2007h . 



Data on the large-scale power spectrum (or correlation function) of galaxies complement the 
CMB measurements by connecting the inhomogeneities observed at redshift z ~ 10 3 in the CMB 
to fluctuations in galaxy density close to z = 0, and by relating fluctuations in gravitating matter 
to fluctuations in luminous matte r (which is an a dditio nal complication). For a recent discussion 
of the galaxy power spectrum see iPercival et al.l (120071 ). from which earlier developments may be 
traced. It is a remarkable success of the current cosmological model that it succeeds in providing a 
reasonable fit to both sets of data. The combination of WMAP data with clustering measurements 
from SDSS or th e 2dFGRS reduces seve ral of the paramet e r unce rtainties. For recent examples of 
such analyses see iTegmark et al.l (|2004al ) and iDoran et al.l (|2007al ) . 



The peak of the galaxy power spectrum reflects the Hubble length at matter-radiation equality 
and so constrains Vty[h. The overall shape of the spectrum is sensitive to the densities of the dif- 
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ferent matter components (e.g., neutrinos would cause damping on small scales) and the density of 
dark energy. The same physics that leads to acoustic peaks in the CMB anisotro py causes oscilla- 
tions in the galaxy power spectrum — or a single peak in the correlation function. Eisenstein et al 



(|2005l ) report a three standard deviation detection of this "baryon acoustic oscillation" peak at 
~ lOO/i -1 Mpc in the correlation function of luminous red galaxies (LRG's) measured in the SDSS. 
The resulting measurement of is independent of and consistent with other low redshift mea- 
surements and with the high redshift WMAP result. This is remarkable given the widely different 
redshifts probed (LRG's probe z = 0.35) and notable because possible systematics are different. 
For discussions of the effi cacy o f futur e measurements of the baryon acous tic oscillation peak t o 
constrain dark energy see IWana (|2006l ) , iMcDonald k, Eisenstein! (|2006l ) , and iDoran et al.l (|2007bl ) . 
For constraints from a join t analysis of these data with supernovae and CMB anisotropy data see 
Wang Mukherieej (hood l. 



Tegmark et al.1 (120061 ) include a nice description of how the large-scale galaxy power spectrum 
provides independent measurement of and Qb, which breaks several parameter degeneracies 
and thereby decreases uncertainties on Om,^ and to- A combined WMAP+SDSS analysis reduces 
uncertainties on the matter d ensity , neutrino density, and tensor-to-scalar ratio by roughly a factor 
of two. See ISanchez et al.l (|2006l ) for an analysis of the 2dFGRS large-scale structure data in 
conjunction with CMB measurements. 

Measurements of the clustering of Lyman-a forest clouds complement larger-scale constraints, 
such as those from the CMB and large-scale st ructure, by probing the power spectrum of fluctu- 



ations on smaller scales (McDonald et al 



200.4) . Combi ning observations of 3000 SDSS Lyman-a 



forest cloud spectra with other data, ISeliak et al.l (J2006J) constrain possible variation with scale of 
the spectral index of the primordial power spectrum and find that Lyman-a cloud clustering may 
indicate a slightly higher power spectrum normalization, erg (the fractional mass density inhomo- 
geneity smoothed over 8/i _1 Mpc), than do the WMAP data alone, or the WMAP data combined 
with large-scale structure measurements. 

The presence of dark energy or non-zero spatial curvature causes time evolution of gravitational 
potentials as CMB photons traverse the Universe from their "emission" at z ~ 10 3 to today. The 
resulting net redshifts or blueshifts of photons cause extra CMB anisotropy, known as the Integrated 
Sachs- Wolfe (ISW) contribution. This contribution has been detected by cross-correlation of CMB 
anisotropy and large-scale struc ture data. The resulting constraints on dark energy are consistent 



with the model discussed above (jBoughn Crittendenll2005l ; iGaztahaga et al.ll2006l . and references 



cited therein). In p rinciple, measure ments of the ISW effect at different redshifts can constrain the 
dark energy model. iPogosianl (J2006J) discusses recent developments and the potential of future ISW 
measurements. 
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7.3. Cosmic Complementarity: Combining Measurements 



The plethora of observational constraints on cos mological parame t ers h as spawned int e rest i n 
stat istical methods for com bining these constraints. iLewis Sz Bridld (|2002l ). I Verde et al.l ((2003), 
and iTegmark et al.l (|2004al ) discuss statistical methods employed in some of the recent analyses 
described above. Use of such advanced statistical techniques is important because of the growing 
number of parameters in current models and possible degenera cies between them i n fitting the 



observ a tional data . Deve l opments may be t r aced back thro ugh lAlam et al. (2007), Zh ang et al 



(120071), IZhao et all (|2007n . Davis et al.1 (|2007n . IWrightJ (|2007j), and lKurek fc Szvdlowskil pOOTn 



To describe large-scale features of the Universe (including CMB anisotropy measured by 
WMAP and some smaller-angular-scale experiments, large-scale structure in the galaxy distribu- 
tion, and the SNIa luminosity-distance relation) the simplest version of the "power-law-spectrum 
spatially- flat ACDM model" requires fitting no fewer than six parameters (jSpergel et al.l 120071 ) : 
nonrelativistic matter density parameter Om, baryon density parameter Ob, Hubble constant Hq, 
amplitude of fluctuations as, optical depth to reionization r, and scalar perturbation index n. This 
model assumes that the primordial fluctuations are Gaussian random phase and adiabatic. As 
suggested by its name, this model further assumes that the primordial fluctuation spectrum is a 
power law (running power-spectral index independent of scale dn/d\nk = 0), the Universe is flat 
(Ok = 0), the bulk of the matter density is CDM (Ocdm = Om — Ob) with no contribution from 
hot dark matter (neutrino density 0„ = 0), and that dark energy in the form of a cosmological 
constant comprises the balance of the mass-energy density (Oa = 1 — Om). Of course, constraints 
on this model assume the validity of the CDM structure formation model. 

Combinations of observations provide improved parameter constraints, typically by breaking 
parameter degeneracies. For example, the constraints from WMAP data alone are relatively weak 
for Hq,Q\ and Ok- Other measurements such as from SNela or galaxy clusters are needed to 
break the degeneracy between Ok and Oa, which lies approximately along Ok ~ —0.3 + 0.4Oa- 
The degeneracy between Om and as is broken by including weak lensing and cluster measurements. 
The degeneracy between Om and Hq can be removed, of course, by including a constraint on Hq. 
As a result, including Hq data restricts the geometry to be very close to flat. A caveat regarding 
this last conclusion is that it assumes that the dark energy density does not evolve. 

CDM-model-dependent clustering limits on baryo n density (Or = (0.02 22 ± 0.0014)/i~ 2 from 
WMAP and SDSS data combined at 95 % confidence, ITegmark et al.l 120061 ) are now better than 
those from light element abundance data (because of the tension between the 4 He and 7 Li data 
and the D data). It is important that the galaxy observations complement the CMB data in such a 
way as to lessen reliance on the assumptions stated above for the "power-law flat ACDM model." 
If the SDSS LRG P(k) measurement is combined with WMAP data, then several of the prior 
assumptions used in the WMAP-alone analysis (Ok = 0,0^ = 0, no running of the spectral index 
n of scalar fluctuations, no inflationary gravity waves, no dark energy temporal evolution) are not 
important. A major reason for this is the sensitivity of the SDSS LRG P{k) to the baryon acoustic 
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scale, which sets a "standard ruler" at low redshift. 



The SNela observations are a powerful complement to CMB anisotropy measurements because 
the degeneracy in versus 0^ for SNIa measurements is almost orthogonal to that of the CMB. 
Without any assumption about the value of the Hubble constant but assuming that the dark energy 
does not evolve, combining SNIa and CMB anisotropy data clearly favors nearly flat cosmologies. 
On the other hand, assuming the Universe is spatially flat, combined SNIa and cluster baryon 



fraction data favors dark ener gy th at does not evo l ve — a cosmological constant — see lRapetti et al 



(|2005lUAlcaniz k Zhul (|2005l b and lWilson et al.1 (j2006h 



The bottom line is that statistical analyses of these complementary observations strongly 
support the flat ACDM cosmological model. It is remarkable that many of the key parameters are 
now known to better than 10%. However, several weaknesses remain, as discussed in the following, 
and final, section of this review. Time and lots of hard work will tell if these weaknesses are simply 
details to be cleaned up, or if they reveal genuine failings of the model, the pursuit of which will 
lead to a deeper understanding of physics and/or astronomy. It is worth recalling that, at the 
beginning of the previous century just before Einstein's burst of 1905 papers, it was thought by 
most physicists that classical physics fit the data pretty well. 



8. Open Questions and Missing Links 



We conclude this review by emphasizing that cosmology is by no means "solved." Here we 
list some outstanding questions, which we do not prioritize, although the first two questions are 
certainly paramount (What is most of the Universe made of?). It may interest the rea der to 
comp are this discussion of outstanding problems in cosmology to those discussed in 1996 (jTurok 
Recent discussions of k ey questions, with regard to fund i ng fo r answering such questions, 



comp a: 
1993) • 



may be found in reports of the lNational Research Council! (|200ll . 120031 ). 



8.1. What is "Dark Energy"? 



As discussed in Sees. [4] and there is strong evidence that the dominant component of mass- 
energy is in the form of something like Einstein's cosmological constant. In detail, does the dark 
energy vary with space or time? Data so far are consistent with a cosmological constant with no 
spatial or temporal evolution, but the constraints do not strongly exclude other possibilities. This 
un certainty is complemented by t he re lativel y weak direct evidence for a spatially-flat universe; 



as 



Wrightl (|2006l ). iTegmark et al.1 (|2006l ). and I Wang &: Mukheried (|2007l ) point out, it is incorrect 
to assume Qk = when constraining the dark energy time dependence, because observational 
evidence for spatial flatness assumes that the dark energy does not evolve. 



More precisely, dark energy is often described by the XCDM parameterization, where it is 
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assumed to be a fluid with pressure px = ^xPx, where px is the energy density and wx is time- 
independent and negative but not ne cessaril y —1 a s in the ACDM model. This is an inaccurate 
parameterization of dark energy; see iRatral (|199ll ) for a discussion of the scalar field case. In 
addition, dark energy and dark matter are coupled in some models now under discussion , so th is 
also needs to be accou n ted for when co mparin g data and mode ls; see lAmendola et al.l (|2007l ). 
Bonometto et al. (j2007l ). Guo et al. ( 2007 ). and Balbi et al. ( 2007 ) for recent discussions. 



On the astronomy side, the evidence is not iron-clad; for example, inference of the presence 
of dark energy from CMB anisotropy data relies on the CDM structure formation model and the 
SNIa redshift-magnitude results require extraordinary nearly "standard candle" -like behavior of 
the objects. Thus, work remains to be done to measure (or reject) dark energy spatial or temporal 
variation and to shore up the observational methods already in use. 

With tighter observational constraints on "dark energy," one might hope to be guided to a 
more fundamental model for this construct. At present, dark energy (as well as dark matter) 
appears to be somewhat disconnected from the rest of physics. 



8.2. What is Dark Matter? 

Astronomical observations currently constrain most of the gravitating matter to be cold (small 
primeval free-streaming velocity) and weakly interacting. Direct detection would be more satisfying 
and this probably falls to laboratory physicists to pursue. The Large Hadron Collider (LHC) may 
produce evidence for the supersymmetric sector that provides some of the most-discussed current 
options for the culprit. As mentioned in the previous question, some models allow for coupling 
between the dark matter and dark energy. On the astronomy side, observations may provide 
further clues and, perhaps already do; there are suggestions of problems with "pure" CDM from 
the properties of dwarf galaxies and galactic nuclear density profiles. Better understanding of the 
complex astrophysics that connect luminous (or, at least, directly detectable) matter to dark matter 
will improve such constraints. 



8.3. What are the Masses of the Neutrinos? 



In contrast to various proposed candidates for the more dominant "cold" component of dark 
matter, we know that ne utrinos exist. Whil e there are indications from underground experiments of 
non-zero neutrino mass (jEguchi et al.ll2003l ) and the cosmological tests discussed above yield upper 
bounds on the sum of masses of all light neutrino species, there has yet to be a detection of the 
effect of ne utrinos on structu re formation. A highly model-dependent analysis of Lyman-a forest 
clustering (jSeljak et al.ll2006l ) results in an upper bound of ^2m u < 0.17 eV (95 % confidence; the 
sum is over light neutrino species). 
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8.4. Are Constraints on Baryon Density Consistent? 



Using the standard theory for nucleosynthesis to constrain the baryon density from observations 
of light element abundances, measureme nts of 4 He and 7 L i impl y a hi gher baryon den sity than do D 
measurements, see Sees. 12.31 and 17.21 and iField h Sarkaii (|2006l ) and ISteigmanl (|2006l ). Constraints 
on the baryon density from WMAP CMB anisotropy data are consistent with that from the D 
abundance measurements. It is possible that more and better data will resolve this discrepancy. 
On the other hand, this might be an indication of new physics beyond the standard model. 



8.5. When and How Was the Baryon Excess Generated? 

Matter is far more common than anti-matter. It is not yet clear how this came to be. One 
much-discussed option is that grand unification at a relatively high temperature is responsible 
for the excess. An alternate possibility is that the matter excess was generated at much lower 
temperature during the electroweak phase transition. 



8.6. What is the Topology of Space? 



The observational constraints we have reviewed are local; they do not constrain the global 
topology of space. On the largest observabl e scales, CMB an isotropy data may be used to constrain 
models for the topology of space (see, e.g., Key et al.1 120071 . and references cited therein). Current 
data do not indicate a real need for going beyond the simplest spatially-flat Euclidean space with 
trivial topology. 



8.7. What Are the Initial Seeds for Structure Formation? 

The exact nature of the primordial fluctuations is still uncertain. The currently-favored expla- 
nation posits an inflationary epoch that precedes the conventional Big Bang era (see Sec. [3]). The 
simplest inflation models produce nearly scale-invariant adiabatic pertur bations. A key const raint 



on inflation models is the slope of the primordial spectrum; WMAP data (ISpergel et al.ll2007l ) sug- 
gest a deviation from the scale-invariant n = 1 value, but this is not yet well measured. At present, 
the most promising method for observationally probing this early epoch is through detection of (the 
scale-invariant spectrum of) inflationary gravity waves predicted in a number of inflation models. 
Detection of these waves or their effects (e.g., measuring the ratio of tensor to scalar fluctuations 
via CMB anisotropy data), would constrain models for inflation; however, non-detection would not 
rule out inflation because there are simple inflation models without significant gravity waves. 



Another critical area for studying the initial fluctuations regards the possibility of non-Gaussian 
perturbations or isocurvature (rather than adiabatic) perturbations. The evidence indicates that 
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these are sub-dominant, but that does not exclude a non-vanishing and interesting contribution. 

Some models of inflation also predict primordial magnetic field fluctuations. These can have 
effects in the low-redshift Universe, including on the CMB anisotropy. Observational detection of 
some of these effects will place interesting constraints on inflation. 



8.8. Did the Early Universe Inflate and Reheat? 

Probably (although we would not be astonished if the answer turned out to be no). With 
tighter observational constraints on the fossil fluctuations generated by quantum mechanics during 
inflation one might hope to be guided to a more fundamental model of inflation. At present, 
inflation is more of a phenomenological construct; an observationally-consistent, more fundamental 
model of inflation could guide the development of very high energy physics. This would be a major 
development. Another pressing need is to have a more precise model for the end of inflation, when 
the Universe reheats and matter and radiation are generated. It is possible that the matter excess 
is generated during this reheating transition. 



8.9. When, How, and What Were the First Structures Formed? 

Discovery of evidence for the epoch of reionization, from observations of absorption line systems 
toward high-redshift quasars and the polarization anisotropy of the CMB, has prompted intense 
interest, both theoretical and observational, in studying formation of the first objects. See Sec. 15.31 
above. 



8.10. How Do Baryons Light Up Galaxies and What Is Their Connection to Mass? 

Carrying on from the previous question, the details of the process of turning this most familiar 
component of mass-energy into stars and related parts of galaxies remains poorly understood. Or 
so it seems when compared with the much easier task of predicting how collisionless dark matter 
clusters in a Universe dominated by dark matter and dark energy. Important problems include the 
effects of "feedback" from star formation and active galactic nuclei, cosmic reionization, radiative 
transfer, and the effect of baryons on halo profiles. High-resolution hydrodynamic simulations are 
getting better, b ut even Moore's law will not help much in the very near future (see comment in 



Gott et al.l 120061 ) . Solving these problems is critical, not only for understanding galaxy formation, 
but also for using galaxies — the "atoms of cosmology" — as a probe of the properties of dark 
matter and dark energy. 



Clues to the relationship between mass and light and, therefore, strong constraints on models 
of galaxy formation, include the detailed dependence of galaxy properties on environment. Out- 
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standing puzzles include the observation that, while galaxy morphology and luminosity strongly 
vary with environment, the properties of early-typ e (elliptical and SO) galaxies (particularly their 
colors) are remarkably insensitive to environment ([Park et al,l 120071 ). 



8.11. How Do Galaxies and Black Holes Coevolve? 



It is now clear that nearly every sufficiently massive galaxy harbors a supermassive black hole 
in its core. The masses of the central supermassive black holes are f ound to correlate strongly 
with properties of th e host galaxy, including bulge velocity dispersion (jFerrarese k, Merrittl 120001 ; 
Gebhardt et al.ll2000l ). Thus, galaxy formation and the formation and feedin g of black holes are 



intim ately related (see, e.g., ISilk &; Rees 
20051 ). 



1998 



Kauffmann &: Haehneltl 12000 : Begelman fe Nath 



8.12. Does the Gaussian, Adiabatic CDM Structure Formation Model Have a Real 

Flaw? 

This model works quite well on large scales. However, on small scales it appears to have too 
much power at low redshift (excessively cuspy halo cores, excessively large galactic central densi- 
ties, and too many low-mass satellites of massive galaxies). Modifications of the power spectrum 
to alleviate this excess small-scale power cause too little power at high redshift and thus delay 
formation of clusters, galaxies, and Lyman-a clouds. Definitive resolution of this issue will require 
more and better observational data as well as improved theoretical modeling. If the CDM structure 
formation model is found to be inadequate, this might have significant implications for a number 
of cosmological tests that assume the validity of this model. 



8.13. Is the Low Quadrupole Moment of the CMB Anisotropy a Problem for Flat 

ACDM? 

The small amplitude of the quadrupole moment observed by COBE persists in the WMAP ob- 



serva tions even after many rounds of reanalysis of possible foreground contributions (see lPark et al 



2007, and references cited therein). Although one cannot, by definition, rule out the possibility 
that it is simply a statistical fluke (with significance of about 95 % in flat ACDM), this anomaly 
inspires searches for alternative models, including multiply-connected Universes (see above). 
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8.14. Are the Largest Observed Structures a Problem For Flat ACDM? 



The largest superclusters, e.g., the "Sloan Great Wall" (jGott et al.l 120051 ). seen in galaxy 
redshift surveys are not reproduced by simulations of the concordance fla t ACDM cosmol ogy 
(jEinasto et al.l 120061 ) . Perhaps we need larger simulations (see discussion in iGott et al.l 120061 ) or 
better understanding of how galaxies trace mass. 



8.15. Why Do We Live Just Now? 

Last, because we see the Universe from only one place, at only one time, we must wrestle with 
questions related to whether or not we (or at least our location) is special. 



Peebles! (|2005l ) notes the remarkable coincidences that we observe the Universe when (1) it 
has just begun making a transition from being dominated by matter to being dominated by dark 
energy, (2) the Milky Way is just running out of gas for forming stars and planetary systems, and 
(3) galaxies have just become useful tracers of mass. While anthropic arguments have been put 
forward to answer the question of why we appear to live at a special time in the history of the 
Universe, a physically motivated answer might be more productive and satisfying. Understanding 
of the details of structure formation, including conversion of baryons to stars (mentioned above), 
and constraints on possible evolution of the components of mass-energy in the Universe may provide 
clues. 

Progress in cosmology is likely to come from more and higher-quality observational and sim- 
ulation data as well as from new ideas. A number of ground-based, space-based, and numerical 
experiments continue to collect data and new near-future particle physics, cosmology, astronomy, 
and numerical experiments are eagerly anticipated. It is less straightforward to predict when a 
significant new idea might emerge. 
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0507463 and the hospitality of the Department of Astrophysical Sciences at Princeton University 
during sabbatical leave. 



-47- 



REFERENCES 

Aarseth, S. J., Turner, E. L., & Gott, III, J. R. 1979, ApJ, 228, 664 
Abbott, L. F., & Wise, M. B. 1984, Nucl. Phys. B, 244, 541 
Abell, G. 0. 1958, ApJS, 3, 211 
Abell, G. 0. 1961, AJ, 66, 607 

Abell, G. O., Corwin, H. G., & Olowin, R. P. 1989, ApJS, 70, 1 
Adelberger, K. L., et al. 2005, ApJ, 619, 697 

Alam, U., Sahni, V., & Starobinsky, A. A. 2007, J. Cosmology Astropart. Phys., 0702, 011 

Albrecht, A., & Steinhardt, P. J. 1982, Phys. Rev. Lett., 48, 1220 

Alcaniz, J. S., Dev, A., & Jain, D. 2005, ApJ, 627, 26 

Alcaniz, J. S., & Zhu, Z.-H. 2005, Phys. Rev. D, 71, 083513 

Allen, S. W., et al. 2004, MNRAS, 353, 457 

Alpher, R. A. 1948, Phys. Rev., 74, 1577 

Alpher, R. A., Follin, Jr., J. W., & Herman, R. 1953, Phys. Rev., 92, 1347 
Alpher, R. A., & Herman, R. 1948, Nature, 162, 774 

Alpher, R. A., & Herman, R. 2001, Genesis of the Big Bang (Oxford: Oxford Univ.) 
Alvarez, M. A., Shapiro, P. R., Ahn, K., & Iliev, I. I. 2006, ApJ, 644, L101 
Amendola, L., Campos, G. C, & Rosenfeld, R. 2007, Phys. Rev. D, 75, 083506 
Arp, H., k Bertola, F. 1969, Astrophys. Lett., 4, 23 
Ashman, K. M. 1992, PASP, 104, 1109 
Astier, P., et al. 2006, A&A, 447, 31 

Avila-Reese, V. 2006, in Solar, Stellar and Galactic Connections between Particle Physics and 
Astrophysics, eds. A. Carramihana, F. S. Guzman, and T. Matos (Berlin: Springer), 115 

Babcock, H. W. 1939, Lick Obs. Bull., 19, 41 

Baccigalupi, C, et al. 2002, Phys. Rev. D, 65, 063520 

Bahcall, N. A., et al. 2003, ApJS, 148, 243 



-48- 



Balbi, A., Bruni, M., & Quercellini, C. 2007, astro-ph/0702423 



Balbi, A., Natoli, P., & Vittorio, N. 2006, in Cosmic Polarization, ed. R. Fabbri (RSP), in press, 
|astro-ph/060651l] 

Baltz, E. A. 2004, eConf C040802, L002 (Available at www.stanford.slac.edu/- 
econf /c040802/proceedings .htm.) 

Bardeen, J. M. 1980, Phys. Rev. D, 22, 1882 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 
Bardeen, J. M., Steinhardt, P. J., & Turner, M. S. 1983, Phys. Rev. D, 28, 679 
Barger, V., Gao, Y., & Marfatia, D. 2007, Phys. Lett. B, 648, 127 
Barnes, J., & Hut, P. 1986, Nature, 324, 446 

Bartolo, N., Komatsu, E., Matarrese, S., & Riotto, A. 2004, Phys. Rep., 402, 103 
Baugh, C. M. 2006, Rept. Prog. Phys., 69, 3101 

Bean, R., Dunkley, J., & Pierpaoli, E. 2006, Phys. Rev. D, 74, 063503 
Becker, R. H., et al. 2001, AJ, 122, 2850 

Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 
Beckwith, S. V. W., et al. 2006, AJ, 132, 1729 
Begelman, M. C, & Nath, B. B. 2005, MNRAS, 361, 1387 
Beichman, C. A., et al. 1988, NASA RP-1190, Vol. 1 

Bekenstein, J. D. 2004, Phys. Rev. D, 70, 083509, erratum 71, 069901 (2005) 
Bennett, C. L., et al. 2003a, ApJS, 148, 1 
Bennett, C. L., et al. 2003b, ApJS, 148, 97 

Berendzen, R., Hart, R., & Seeley, D. 1976, Man Discovers the Galaxies (New York: Science 
History) 

Berger, M. J., & Collela, B. 1989, J. Comp. Phys., 82, 64 
Berlind, A. A., & Weinberg, D. H. 2002, ApJ, 575, 587 

Bernardeau, F., Colombi, S., Gaztahaga, E., & Scoccimarro, R. 2002, Phys. Rep., 367, 1 
Bernstein, J., Brown, L. S., & Feinberg, G. 1989, Rev. Mod. Phys., 61, 25 



- 49 - 



Bernstein, J., k Feinberg, G., eds. 1986, Cosmological Constants: Papers in Modern Cosmology 
(New York: Columbia Univ.) 

Bertone, C, Hooper, D., k Silk, J. 2005, Phys. Rep., 405, 279 

Bertotti, B., Balbinot, R., Bergia, S., k Messina, A., eds. 1990, Modern Cosmology in Retrospect 
(Cambridge: Cambridge Univ.) 

Bertschinger, E. 1998, ARA&A, 36, 599 

Bertschinger, E., k Dekel, A. 1989, ApJ, 336, L5 

Bertschinger, E., k Gelb, J. M. 1991, Comp. Phys., 5, 164, 178, 179 

Binney, J. 1977, ApJ, 215, 483 

Blumenthal, G. R., Faber, S. M., Primack, J. R., k Rees, M. J. 1984, Nature, 311, 517 
Blumenthal, G. R., Pagels, H., k Primack, J. R. 1982, Nature, 299, 37 
Bode, P., Ostriker, J. P., k Xu, G. 2000, ApJS, 128, 561 
Bohringer, H., et al. 2004, A&A, 425, 367 

Bond, J. R., Efstathiou, G., k Silk, J. 1980, Phys. Rev. Lett., 45, 1980 
Bond, J. R., Szalay, A. S., k Turner, M. S. 1982, Phys. Rev. Lett., 48, 1636 
Bonnor, W. B. 1957, MNRAS, 117, 104 

Bonometto, S. A., Casarini, L., Colombo, L. P. L., k Mainini, R. 2007, |astro-ph/06 12672"] 

Borgani, S. 2006, in Lecture Notes in Physics, 2005 Guillermo Haro Summer School on Clusters 
(Springer), in press, [astro-ph /0605575 

Boughn, S. P., k Crittenden, R. G. 2005, New A Rev., 49, 75 

Branch, D. 1998, ARA&A, 36, 17 

Brax, P., Martin, J., k Riazuelo, A. 2000, Phys. Rev. D, 62, 103505 

Broadhurst, T. J., Ellis, R. S., Koo, D. C, k Szalay, A. S. 1990, Nature, 343, 726 

Brodie, J. P., k Strader, J. 2006, ARA&A, 44, 193 

Bromm, V., k Larson, R. B. 2004, ARA&A, 42, 79 

Brout, R., Englert, F., k Gunzig, E. 1978, Ann. Phys., 115, 78 

Brown, I., k Crittenden, R. 2005, Phys. Rev. D, 72, 063002 



-50- 



Buonanno, A. 2004, in Particle Physics and Cosmology: The Quest for Physics Beyond the Standard 
Model(s), eds. H. E. Haber and E. Nelson (Singapore: World Scientific), 855 

Caldwell, R. R., & Doran, M. 2004, Phys. Rev. D, 69, 103517 

Carlstrom, J. E., Holder, G. P., & Reese, E. D. 2002, ARA&A, 40, 643 

Carr, B. 1994, ARA&A, 32, 531 

Carroll, S. M. 2004, in Measuring and Modeling the Universe, ed. W. L. Freeman (Cambridge: 
Cambridge Univ.), 235 

Cen, R. 1992, ApJS, 78, 341 

Chae, K.-H., et al. 2002, Phys. Rev. Lett., 89, 151301 

Chae, K.-H., Chen, C, Ratra, B., k Lee, D.-W. 2004, ApJ, 607, L71 

Challinor, A. 2005, in The Physics of the Early Universe, ed. E. Papantonopoulos (Berlin: Springer), 
71 



Chandrasekhar, S., & Henrich, L. R. 1942, ApJ, 95, 
Chang, Z., Wu, F.-Q., & Zhang, X. 2006, Phys. Lett. B, 633, 14 
Charlier, C. V. L. 1925, PASP, 37, 115 

Chen, C, Gott, III, J. R., & Ratra, B. 2003, PASP, 115, 1269 
Chen, C, & Ratra, B. 2003a, ApJ, 582, 586 
Chen, C, & Ratra, B. 2003b, PASP, 115, 1143 
Chen, C, & Ratra, B. 2004, ApJ, 612, LI 

Choudhury, T. R., & Ferrara, A. 2006a, in Cosmic Polarization, ed. R. Fabbri (RSP), in press, 
|astro-ph/0603149] 

Choudhury, T. R., & Ferrara, A. 2006b, MNRAS, 371, L55 

Cline, J. M. 2006, in Les Houches Summer School, Session 86: Particle Physics and Cosmology: 
The Fabric of Spacetime, hep-ph/0609145| 

Clocchiatti, A., et al. 2006, ApJ, 642, 1 

Clowe, D., et al. 2006, ApJ, 648, L109 

Coil, A. L., et al. 2004, ApJ, 609, 525 

Cole, S., et al. 1994, MNRAS, 271, 781 



-51- 



Coleman, S., & De Luccia, F. 1980, Phys. Rev. D, 21, 3305 
Colless, M., et al. 2001, MNRAS, 328, 1039 
Constantin, A., & Vogeley, M. S. 2006, ApJ, 650, 727 
Cooray, A., & Sheth, R. K. 2002, Phys. Rep., 372, 1 

Copeland, E. J., Sami, M., & Tsujikawa, S. 2006, Int. J. Mod. Phys. D, 15, 1753 
Couchman, H. M. P. 1991, ApJ, 368, L23 

Courteau, S., & Willick, J., eds. 2000, Cosmic Flows Workshop (San Francisco: Astron. Soc. Pac.) 
Croom, S. M., et al. 2004, MNRAS, 349, 1397 
Cyburt, R. H. 2004, Phys. Rev. D, 70, 023505 
da Angela, J., et al. 2005, MNRAS, 360, 1040 
da Costa, L. N., et al., 1998, AJ, 116, 1 

Dalton, G. B., Efstathiou, C, Maddox, S. J., & Sutherland, W. J. 1992, ApJ, 390, LI 
Daly, R. A., & Djorgovski, S. G. 2006, |astro-ph/0609791| 

Davis, M., Efstathiou, C, Frenk, C. S., & White, S. D. M. 1985, ApJ, 292, 371 
Davis, M., & Geller, M. J. 1976, ApJ, 208, 13 
Davis, M., k Peebles, P. J. E. 1983, ARA&A, 21, 109 

Davis, T. M., & Lineweaver, C. H. 2004, Pub. Astron. Soc. Australia, 21, 97 
Davis, T. M., et al. 2007, ApJ, 666, 716 

DeGrasse, R. W., Hogg, D. C, Ohm, E. A., & Scovil, H. E. D. 1959, Proc. Natl. Electronics Conf., 
15, 370 

Dekel, A. 1994, ARA&A, 32, 371 

Dekel, A., & Ostriker, J. P., eds. 1999, Formation of Structure in the Universe (Cambridge: Cam- 
bridge Univ.) 

de Lapparent, V., Geller, M. J., & Huchra, J. P. 1986, ApJ, 302, LI 
de Sitter, W. 1917, MNRAS, 78, 3 
de Vaucouleurs, G. 1953, AJ, 58, 30 



-52- 



de Vaucouleurs, G. 1969, Astrophys. Lett., 4, 17 

Diaz-Rivera, L. M., Samushia, L., k Ratra, B. 2006, Phys. Rev. D, 73, 083503 
Dicke, R. H. 1970, Gravitation and the Universe (Philadelphia: American Philo. Soc.) 
Dicke, R. H., k Peebles, P. J. 1965, Space Sci. Rev., 4, 419 

Dicke, R. H., k Peebles, P. J. E. 1979, in General Relativity: An Einstein Centenary Survey, eds. 
S. W. Hawking and W. Israel (Cambridge: Cambridge Univ.), 504 

Dicke, R. H., Peebles, P. J. E., Roll, P. G., k Wilkinson, D. T. 1965, ApJ, 142, 414 

Dine, M., k Kusenko, A. 2004, Rev. Mod. Phys., 76, 1 

Dolgov, A. D. 2002, Phys. Rep., 370, 333 

Doran, M., Robbers, G., k Wetterich, C. 2007a, Phys. Rev. D, 75, 023003 
Doran, M., Stern, S., k Thommes, E. 2007b, J. Cosmology Astropart. Phys., 0704, 015 
Doroshkevich, A. G., k Novikov, I. D. 1964, English translation in 1964, Sov. Phys.-Dokl., 9, 111 
Doroshkevich, A. G., k Novikov, I. D. 1970, English translation in 1971, Sov. Astron.-AJ, 14, 763 
Dressier, A., et al. 1987, ApJ, 313, L37 

Dubinski, J., Kim, J., Park, C, k Humble, R. 2004, New A, 9, 111 
Durrer, R., Novosyadlyj, B., k Apunevych, S. 2003, ApJ, 583, 33 

Eddington, A. S. 1923, The Mathematical Theory of Relativity (Cambridge: Cambridge Univ.), 
161-2 

Efstathiou, G. 1990, in Physics of the Early Universe, eds. J. A. Peacock, A. F. Heavens, and A. 
T. Davies (Adam Hilger, Bristol), 361 

Eguchi, K., et al. 2003, Phys. Rev. Lett., 90, 021802 

Einasto, J. 2005, in Frontiers of the Universe: Cosmology 2003, eds. A. Blanchard and M. Signore 
(Berlin: Springer), 241 

Einasto, J., et al. 2006, A&A, 459, LI 

Einasto, J., Kaasik, A., k Saar, E. 1974, Nature, 250, 309 

Einstein, A. 1917, English translation in The Principle of Relativity, eds. H. A. Lorentz, A. Einstein, 
H. Minkowski, and H. Weyl (Dover, New York, 1952), 175 

Einstein, A., k de Sitter, W. 1932, Proc. Natl. Acad. Sci. (USA), 18, 213 



-53- 



Eisenstein, D. J., et al. 2005, ApJ, 633, 560 

El-Ad, H., Piran, T., & da Costa, L. N. 1996, ApJ, 462, L13 

Elgar0y, 0., & Lahav, O. 2005, New J. Phys., 7, 61 

Ellis, G. F. R., & Stoeger, W. 1988, Class. Quantum Grav., 5, 207 

Ellis, R. S. 2007, |astro-ph/0701024| 

Esmailzadeh, R., Starkman, G. D., &; Dimopoulos, S. 1991, ApJ, 378, 504 

Evrard, A. E. 1988, MNRAS, 235, 911 

Faber, S. M., & Gallagher, J. S. 1979, ARA&A, 17, 135 

Falco, E. E., et al. 1999, PASP, 111, 438 

Fan, X., Carilli, C. L., & Keating, B. 2006, ARA&A, 44, 415 

Fan, X., et al. 2002, AJ, 123, 1247 

Fedeli, C, & Bartelmann, M. 2007, A&A, 461, 49 

Ferguson, H. C, Dickinson, M., & Williams, R. 2000, ARA&A, 38, 667 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 

Field, B. D., & Sarkar, S. 2006, J. Phys. G, 33, 220 (The latest version is available at 
|http://pdg.lbl. gov/1 ) 

Finzi, A. 1963, MNRAS, 127, 21 

Firmani, C, et al. 2001, MNRAS, 321, 713 

Fischler, W., Ratra, B., & Susskind, L. 1985, Nucl. Phys. B, 259, 730 
Fisher, K. B., et al. 1993, ApJ, 402, 42 
Fisher, K. B., et al. 1995, ApJS, 100, 69 
Fixsen, D. J., et al. 1996, ApJ, 473, 576 
Freedman, W. L., et al. 2001, ApJ, 553, 47 

Friedmann, A. 1922, English translation in Bernstein, J., & Feinberg, G., eds. 1986, 49 
Friedmann, A. 1924, English translation in Bernstein, J., & Feinberg, G., eds. 1986, 59 
Fukugita, M., Futamase, T., & Kasai, M. 1990, MNRAS, 246, 24P 



-54- 



Fukugita, M., Hogan, C. J., k Peebles, P. J. E. 1998, ApJ, 503, 518 

Gal, R. R., et al. 2004, AJ, 128, 3082 

Gamow, G. 1946, Phys. Rev., 70, 572 

Gamow, G. 1948, Phys. Rev., 74, 505 

Gamow, G., k Teller, E. 1939, Phys. Rev., 55, 654 

Gawiser, E. 2006, in New Horizons in Astronomy: Frank N. Bash Symposium 2005, eds. S. Kan- 
nappan, et al. (San Francisco: Astron. Soc. Pac), 177 

Gaztahaga, E., Manera, M., k Multamaki, T. 2006, MNRAS, 365, 171 

Gebhardt, K., et al. 2000, ApJ, 539, L13 

Geller, M. J., et al. 1997, AJ, 114, 2205, 2836 

Giovanelli, R., k Haynes, M. P. 1985, AJ, 90, 2445 

Giovanelli, R., k Haynes, M. P. 1991, ARA&A, 29, 499 

Giovannini, M. 2004, Int. J. Mod. Phys. D, 13, 391 

Giovannini, M. 2005, Int. J. Mod. Phys. D, 14, 363 

Giovannini, M. 2006, Class. Quantum Grav., 23, R1-R44 

Gladders, M. D., k Yee, H. K. C. 2005, ApJS, 157, 1 

Gnedin, N. Y., k Bertschinger, E. 1996, ApJ, 470, 115 

Gnedin, N. Y., k Fan, X. 2006, ApJ, 648, 1 

Goldberg, D. M., k Spergel, D. N. 2000, ApJ, 544, 21 

Gorski, K. M., et al. 1998, ApJS, 114, 1 

Gott, III, J. R. 1977, ARA&A, 15, 235 

Gott, III, J. R. 1982, Nature, 295, 304 

Gott, III, J. R., et al. 2006, ApJ, submitted, [astro-ph/0610762 
Gott, III, J. R., et al. 2005, ApJ, 624, 463 
Gott, III, J. R., et al. 1989, ApJ, 340, 625 

Gott, III, J. R., Vogeley, M. S., Podariu, S., k Ratra, B. 2001, ApJ, 549, 1 



-55- 



Gott, III, J. R., Weinberg, D. H., & Melott, A. L. 1987, ApJ, 319, 1 

Gratton, R. G., et al. 2003, A&A, 408, 529 

Gregory, P. C, & Condon, J. J. 1991, ApJS, 75, 1011 

Gregory, S. A., & Thompson, L. A. 1978, ApJ, 222, 784 

Groth, E. J., & Peebles, P. J. E. 1977, ApJ, 217, 385 

Gunn, J. E., & Gott, III, J. R. 1972, ApJ, 176, 1 

Gunn, J. E., k, Peterson, B. A. 1965, ApJ, 142, 1633 

Gunn, J. E., et al. 2006, AJ, 131, 2332 

Guo, Z.-K., Ohta, N., & Tsujikawa, S. 2007, Phys. Rev. D, 76, 023508 

Gush, H. P., Halpern, M., & Wishnow, E. H. 1990, Phys. Rev. Lett., 65, 537 

Guth, A. H. 1981, Phys. Rev. D, 23, 347 

Guth, A. H., & Pi, S.-Y. 1982, Phys. Rev. Lett., 49, 1110 

Guyot, M., & Zel'dovich, Ya. B. 1970, A&A, 9, 227 

Guzzo, L., et al. 1997, ApJ, 489, 37 

Hamilton, A. J. S. 1988, ApJ, 331, L59 

Hamilton, A. J. S. 1998, in The Evolving Universe, ed. D. Hamilton (Dordrecht: Kluwer Academic), 
185 

Hannestad, S. 2006, Annu. Rev. Nucl. Part. Sci., 56, 137 
Hansen, B. M. S., et al. 2004, ApJS, 155, 551 
Harrison, E. R. 1968, Phys. Rev., 167, 1170 
Harrison, E. R. 1970, Phys. Rev. D, 1, 2726 
Harrison, E. 1993, ApJ, 403, 28 

Harrison, E. 2001, Cosmology: The Science of the Universe, 2nd ed. (corrected) (Cambridge: 
Cambridge Univ.) 

Hartwick, F. D. A., k Schade, D. 1990, ARA&A, 28, 437 

Hatton, S. 1999, MNRAS, 310, 1128 



-56- 



Hauser, M. G., & Peebles, P. J. E. 1973, ApJ, 185, 757 
Hawking, S. W. 1982, Phys. Lett. B, 115, 295 
Hawkins, E., et al. 2003, MNRAS, 346, 78 
Hayashi, C. 1950, Prog. Theor. Phys., 5, 224 
Hernquist, L., & Katz, N. 1989, ApJS, 70, 419 
Hetterscheidt, M., et al. 2007, A&A, 468, 859 
Hinshaw, G., et al. 2007, ApJS, 170, 288 

Hockney, R. W., & Eastwood, J. W. 1988, Computer Simulation Using Particles (Bristol: Adam 
Hilger) 



Hogg, D. W. 1999, astro-ph/9905116 



Hoyle, F. 1953, ApJ, 118, 513 

Hoyle, F., & Narlikar, J. V. 1962, Proc. R. Soc. London A, 270, 334 

Hoyle, F., & Tayler, R. J. 1964, Nature, 203, 1108 

Hoyle, F., & Vogeley, M. S. 2002, ApJ, 566, 641 

Hoyle, F., & Vogeley, M. S. 2004, ApJ, 607, 751 

Hu, W., & Dodelson, S. 2002, ARA&A, 40, 171 

Hubble, E. P. 1925, PASA, 5, 261 

Hubble, E. 1926, ApJ, 64, 321. 

Hubble, E. 1929, Proc. Natl. Acad. Sci. (USA), 15, 168 

Huchra, J. 2007, a list of measurements of Hq available at 



http : //cf a-www . harvard . edu/ ~rmchra/| 



Huchra, J. P., Davis, M., Latham, D., & Tonry, J. 1983, ApJS, 52, 89 

Hudson, M. J., Colless, M., Dressier, A., &; Giovanelli, R. 2000, in Courteau, S., & Willick, J., eds. 
2000, 159 

Hudson, M. J., Smith, R. J., Lucey, J. R., & Branchini, E. 2004, MNRAS, 352, 61 
Humason, M. L., Mayall, N. U., & Sandage, A. R. 1956, AJ, 61, 97 
Imbriani, G., et al. 2004, A&A, 420, 625 



-57- 



Jacoby, G. H., et al. 1992, PASP, 104, 599 
Jarrett, T. H., et al. 2000, AJ, 119, 2498 



Jassal, H. K., Bagla, J. S., & Padmanabhan, T. 2006, astro-ph/0601389 



Jeans, J. 1902, Philo. Trans. R. Soc. London A, 199, 1 

Jeans, J. H. 1923, MNRAS, 84, 60 

Jimenez, R., & Loeb, A. 2002, ApJ, 573, 37 

Joeveer, M., Einasto, J., & Tago, E. 1978, MNRAS, 185, 357 

Johnson, G. 2005, Miss Leavitt's Stars: The Untold Story of the Woman Who Discovered How to 
Measure the Universe (New York: Norton) 

Jones, B. J. T., Martinez, V. J., Saar, E., & Trimble, V. 2005, Rev. Mod. Phys., 76, 1211 

Jones, D. H., et al. 2004, MNRAS, 355, 747 

Kahn, F. D., k Woltjer, L. 1959, ApJ, 130, 705 

Kahniashvili, T., & Ratra, B. 2005, Phys. Rev. D, 71, 103006 

Kahniashvili, T., & Ratra, B. 2007, Phys. Rev. D, 75, 023002 

Kahniashvili, T., von Toerne, E., Arhipova, N. & Ratra, B. 2005 Phys. Rev. D, 71, 125009 
Kaiser, N. 1984, ApJ, 284, L9 
Kaiser, N. 1987, MNRAS, 227, 1 

Kamionkowski, M., & Kosowsky, A. 1999, Annu. Rev. Nucl. Part. Sci., 49, 77 

Kamionkowski, M., Ratra, B., Spergel, D. N., & Sugiyama, N. 1994b, ApJ, 434, LI 

Kamionkowski, M., Spergel, D. N., & Sugiyama, N. 1994a, ApJ, 426, L57 

Katgert, P., et al. 1996, A&A, 310, 8 

Kauffmann, G., & Haehnelt, M. 2000, MNRAS, 311, 576 

Kauffmann, G., White, S. D. M., & Guiderdoni, B. 1993, MNRAS, 264, 201 

Kazanas, D. 1980, ApJ, 241, L59 

Key, J. S., Cornish, N. J., Spergel, D. N., & Starkman, G. D. 2007, Phys. Rev. D, 75, 084034 
Kirshner, R. P., Oemler, Jr., A., Schechter, P. L., & Shectman, S. A. 1981, ApJ, 248, L57 



-58- 



Kitching, T. D., et al. 2007, MNRAS, 376, 771 

Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, 522, 82 

Kochanek, C. S. 2006, in Kochanek, C. S., Schneider, P., & Wambsganss, J. 2006, 91 

Kochanek, C. S., Schneider, P., & Wambsganss, J. 2006, Gravitational Lensing: Strong, Weak 
and Micro: Saas-Fee Advanced Course 33, eds. G. Meylan, P. Jetzer, and P. North (Berlin: 
Springer) 

Kolb, E. W., & Turner, M. S. 1990, The Early Universe (Redwood City: Addison- Wesley) 
Komatsu, E., et al. 2003, ApJS, 148, 119 
Kovac, J. M., et al. 2002, Nature, 420, 772 

Kragh, H. 1996, Cosmology and Controversy: The Historical Development of Two Theories of the 
Universe (Princeton: Princeton Univ.) 

Krauss, L. M., & Chaboyer, B. 2003, Science, 299, 65 

Kravtsov, A. V., Klypin, A. A., & Khokhlov, A. M. 1997, ApJS, 111, 73 

Kravtsov, A. V., et al. 2004, ApJ, 609, 35 

Kravtsov, A. V., Nagai, D., & Vikhlinin, A. A. 2005, ApJ, 625, 588 
Kurek, A., & Szydlowski, M. 2007, |astro-ph/0702484| 
Lacey, C, & Cole, S. 1993, MNRAS, 262, 627 

Lahav, O., & Liddle, A. R. 2006, J. Phys. G, 33, 224 (The latest version is available at 

|http7// P dg.ibi.gov7[ ) 

Lahav, O., & Suto, Y. 2004, Living Rev. Relativity, 7, 8 
Lauer, T. R., & Postman, M. 1994, ApJ, 425, 418 
Leavitt, H. S. 1912, Harvard Coll. Obs. Circ, 173, 1 
Le Fevre, O., et al. 2005, A&A, 439, 877 
Leibundgut, B. 2001, ARA&A, 39, 67 

Lemaitre, G. 1925, J. Math. Phys. (Cambridge, MA), 4, 188 

Lemaitre, G. 1927, abridged and updated English translation in 1931, MNRAS, 91, 483 
Lesgourgues, J., & Pastor, S. 2006, Phys. Rep., 429, 307 



-59- 



Lewis, A. 2004, Phys. Rev. D, 70, 043011 

Lewis, A., & Bridle, S. 2002, Phys. Rev. D, 66, 103511 

Liddle, A. R., & Lyth, D. H. 2000, Cosmological Inflation and Large-Scale Structure (Cambridge: 
Cambridge Univ.) 

Lifshitz, E. 1946, English translation in 1946, J. Phys., 10, 116 

Limber, D. N. 1954, ApJ, 119, 655 

Linde, A. D. 1982, Phys. Lett. B, 108, 389 

Lineweaver, C. H., & Davis, T. M. 2005, Sci. Am., 292 (3), 36 

Loeb, A. 2006a, |astro-ph/0603360| 

Loeb, A. 2006b, Sci. Am., 295 (5), 46 

Longair, M. S. 2006, The Cosmic Century: A History of Astrophysics and Cosmology (Cambridge: 
Cambridge) 

Longair, M. S. 2008, Galaxy Formation (Berlin: Springer) 

Loveday, J., Peterson, B. A., Maddox, S. J., & Efstathiou, G. 1996, ApJS, 107, 201 
Lucchin, F., & Matarrese, S. 1985, Phys. Rev. D, 32, 1316 

Lumsden, S. L., Nichol, R. C, Collins, C. A., & Guzzo, L. 1992, MNRAS, 258, 1 
Lynds, R. 1971, ApJ, 164, L73, plate L3 

Maddox, S. J., Efstathiou, G., Sutherland, W. J., & Loveday, J. 1990, MNRAS, 242, 43P 
Maggiore, M. 2000, Phys. Rep., 331, 283 

Mainini, R., Maccio, A. V., Bonometto, S. A., & Klypin, A. 2003, ApJ, 599, 24 
Maio, U., et al. 2006, MNRAS, 373, 869 

Mandelbrot, B. B. 1982, The Fractal Geometry of Nature (San Francisco: Freeman) 
Martin, D. C, et al. 2005, ApJ, 619, LI 

Martinez, V. J., Lopez-Marti, B., & Pons-Borderfa, M.-J. 2001, ApJ, 554, L5 

Martinez, V. J., & Saar, E. 2002, Statistics of the Galaxy Distribution (Boca Raton: Chapman & 
Hall/CRC) 

Massey, R., et al. 2007, Nature, 445, 286 



- 60 - 



Mather, J. C, et al. 1999, ApJ, 512, 511 

Maurogordato, S., & Lachieze-Rey, M. 1987, ApJ, 320, 13 

McCrea, W. H. 1968, Science, 160, 1295 



McDonald, P., & Eisenstein, D. J. 2006, Phys. Rev. D, submitted, [astro-ph/ 0607122 
McDonald, P., et al. 2005, ApJ, 635, 761 
McKellar, A. 1941, Publ. Dominion Astrophys. Obs., 7, 251 
Merz, H., Pen, U.-L., & Trac, H. 2005, New A, 10, 393 
Michie, R. W. 1969, Kitt Peak Natl. Obs. Contr., 440, 1 
Milgrom, M. 1983, ApJ, 270, 365 
Milgrom, M. 2002, Sci. Am., 287 (2), 42 
Milne, E. A. 1933, ZAp, 6, 1 
Misner, C. W. 1969, Phys. Rev. Lett., 22, 1071 

Misner, C. W., Thorne, K. S., & Wheeler, J. A. 1973, Gravitation (San Francisco: Freeman) 

Mo, H. J., & White, S. D. M. 1996, MNRAS, 282, 347 

Monaghan, J. J. 1992, ARA&A, 30, 543 

Moore, B., et al. 1999a, ApJ, 524, L19 

Moore, B., et al. 1999b, MNRAS, 310, 1147 

Mota, D. F., & van de Bruck, C. 2004, A&A, 421, 71 

Mowbray, A. G. 1938, PASP, 50, 275 

Mukhanov, V. 2005, Physical Foundations of Cosmology (Cambridge: Cambridge Univ.) 

Mukhanov, V. F., Feldman, H. A., & Brandenberger, R. H. 1992, Phys. Rep., 215, 203 

Mukherjee, P., et al. 2003b, ApJ, 598, 767 

Mukherjee, P., et al. 2003a, ApJ, 592, 692 

Mukherjee, P., et al. 2003c, Int. J. Mod. Phys. A, 18, 4933 

Muhoz, C. 2004, Int. J. Mod. Phys. A, 19, 3093 



Munshi, D., Valageas, P., van Waerbeke, L., & Heavens, A. 2006, astro-ph/0612667 



-61- 



Narlikar, J. V., & Padmanabhan, T. 1991, ARA&A, 29, 325 
Narlikar, J. V., et al. 2003, ApJ, 585, 1 

National Research Council 2001, Astronomy and Astrophysics in the New Millenium (Washington: 
National Academy) 

National Research Council 2003, Connecting Quarks with the Cosmos: Eleven Science Questions 
for the New Century (Washington: National Academy) 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Nesseris, S., & Perivolaropoulos, L. 2005, Phys. Rev. D, 72, 123519 



Nobbenhuis, S. 2006, Utrecht Univ. Ph.D. thesis, gr-qc/0609011 
Norberg, P., et al. 2002, MNRAS, 332, 827 

Norman, M. L., & Bryan, G. L. 1999, in Numerical Astrophysics, eds. S. M. Miyama, K. Tomisaka, 
and T. Hanawa (Boston: Kluwer), 19 

North, J. D. 1990, The Measure of the Universe: A History of Modern Cosmology (New York: 
Dover) 

Novikov, I. D. 1964, English translation in 1964, Sov. Phys. JETP, 19, 467 
Ohm, E. A. 1961, Bell Syst. Tech. J., 40, 1065 
Oort, J. H. 1940, ApJ, 91, 273 
Opik, E. 1922, ApJ, 55, 406 
Osmer, P. S. 1981, ApJ, 247, 762 

Ostriker, J. P., & Peebles, P. J. E. 1973, ApJ, 186, 467 
Ostriker, J. P., Peebles, P. J. E., & Yahil, A. 1974, ApJ, 193, LI 
Outram, P. J., et al. 2003, MNRAS, 342, 483 

Padmanabhan, T. 1993, Structure Formation in the Universe (Cambridge: Cambridge Univ.) 
Padmanabhan, T. 2005, Curr. ScL, 88, 1057 
Page, L. A., et al. 2007, ApJS, 170, 335 
Page, L. A., et al. 2003, ApJS, 148, 233 
Page, T. 1952, ApJ, 116, 63 



-62- 



Park, C, et al. 2007, ApJ, 658, 898 

Park, C, Vogeley, M. S., Geller, M. J., k Huchra, J. P. 1994, ApJ, 431, 569 

Park, C.-G., Park, C, k Gott, III, J. R. 2007, ApJ, 660, 959 

Park, C.-G., Park, C., k Ratra, B. 2002, ApJ, 568, 9 

Park, C.-G., Park, C., Ratra, B., k Tegmark, M. 2001, ApJ, 556, 582 

Partridge, R. B. 1995, 3 K: The Cosmic Microwave Background Radiation (Cambridge: Cambridge 
Univ.) 

Peebles, P. J. E. 1965, ApJ, 142, 1317 

Peebles, P. J. E. 1966a, Phys. Rev. Lett., 16, 410 

Peebles, P. J. E. 1966b, ApJ, 146, 542 

Peebles, P. J. E. 1967, Primeval Galaxies, unpublished manuscript based on a paper presented at 
the 4 th Texas Conference on Relativistic Astrophysics (May 1967) 

Peebles, P. J. E. 1968, ApJ, 153, 1 

Peebles, P. J. E. 1971, Physical Cosmology (Princeton: Princeton Univ.) 

Peebles, P. J. E. 1980, The Large-Scale Structure of the Universe (Princeton: Princeton Univ.) 
Peebles, P. J. E. 1982, ApJ, 263, LI 
Peebles, P. J. E. 1984, ApJ, 284, 439 
Peebles, P. J. E. 1990, ApJ, 362, 1 

Peebles, P. J. E. 1993, Principles of Physical Cosmology (Princeton: Princeton Univ.) 

Peebles, P. J. E. 1999, in Dekel, A., k Ostriker, J. P., eds. 1999, 435 

Peebles, P. J. E. 2005, Nucl. Phys. Proc. Suppl., 138, 5 

Peebles, P. J. E., k Dicke, R. H. 1968, ApJ, 154, 891 

Peebles, P. J. E., k Ratra, B. 1988, ApJ, 325, L17 

Peebles, P. J. E., k Ratra, B. 2003, Rev. Mod. Phys., 75, 559 

Peebles, P. J. E., k Yu, J. T. 1970, ApJ, 162, 815 

Pen, U.-L. 1997, New A, 2, 309 



-63- 



Penzias, A. A., & Wilson, R. W. 1965, ApJ, 142, 419 



Percival, W. J. 2006, astro-ph/0601538 



Percival, W. J., et al. 2007, ApJ, 657, 645 

Perivolaropoulos, L. 2006, in Recent Advances in Astronomy and Astrophysics, ed. N. Solomos 
(Melville: AIP), 698 

Perlmutter, S., et al. 1999, ApJ, 517, 565 



Perlmutter, S., et al. 2006, at http://snap.lbl.gov/ 



Pike, R. W., & Hudson, M. J. 2005, ApJ, 635, 11 

Podariu, S., Daly, R. A., Mory, M. P., & Ratra, B. 2003, ApJ, 584, 577 

Podariu, S., Nugent, P., & Ratra, B. 2001, ApJ, 553, 39 

Podariu, S., & Ratra, B. 2000, ApJ, 532, 109 

Podariu, S., et al. 2001, ApJ, 559, 9 

Pogosian, L. 2006, New A Rev., 50, 932 

Pope, A. C, et al. 2004, ApJ, 607, 655 

Postman, M., Huchra, J. P., & Geller, M. J. 1992, ApJ, 384, 404 

Postman, M., et al. 1996, AJ, 111, 615, 988 

Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 

Primack, J. R. 2005, New A Rev., 49, 25 

Quinn, H. R., & Witherell, M. S. 1998, Sci. Am., 279 (4), 76 

Rapetti, D., Allen, S. W., & Weller, J. 2005, MNRAS, 360, 555 

Ratcliffe, A., et al. 1998, MNRAS, 300, 417 

Ratra, B. 1985, Phys. Rev. D, 31, 1931 

Ratra, B. 1988, Phys. Rev. D, 38, 2399 

Ratra, B. 1991, Phys. Rev. D, 43, 3802 

Ratra, B. 1992a, ApJ, 391, LI 

Ratra, B. 1992b, Phys. Rev. D, 45, 1913 



-64- 



Ratra, B., & Peebles, P. J. E. 1988, Phys. Rev. D, 37, 3406 

Ratra, B., & Peebles, P. J. E. 1994, ApJ, 432, L5 

Ratra, B., & Peebles, P. J. E. 1995, Phys. Rev. D, 52, 1837 

Ratra, B., & Quillen, A. 1992, MNRAS, 259, 738 

Rauch, M. 1998, ARA&A, 36, 267 

Rees, M. J., & Ostriker, J. P. 1977, MNRAS, 179, 541 

Refregier, A. 2003, ARA&A, 41, 645 

Refregier, A., et al. 2006, Proc. SPIE, 6265, 62651Y 

Riess, A. G., et al. 1998, AJ, 116, 1009 

Riess, A. G., et al. 2007, ApJ, 659, 98 

Rindler, W. 1956, MNRAS, 116, 662 

Robertson, H. P. 1928, Philo. Mag., 5, 835 

Robertson, H. P. 1929, Proc. Natl. Acad. Sci. (USA), 15, 822 

Romer, A. K., et al. 1994, Nature, 372, 75 

Rosati, P., Borgani, S., & Norman, C. 2002, ARA&A, 40, 539 

Rowan-Robinson, M. 1985, The Cosmological Distance Ladder: Distance and Time in the Universe 
(New York: Freeman) 

Rubakov, V. A., Sazhin, M. V., & Veryaskin, A. V. 1982, Phys. Lett. B, 115, 189 

Rubin, V. C. 1954, Proc. Natl. Acad. Sci. (USA), 40, 541 

Rubin, V. C, & Ford, Jr., W. K. 1970, ApJ, 159, 379 

Rubin, V. C, Thonnard, N., Ford, Jr., W. K., & Roberts, M. S. 1976, AJ, 81, 719 

Ryu, D., Ostriker, J. P., Kang, H., & Cen, R. 1993, ApJ, 414, 1 

Sahni, V., & Coles, P. 1995, Phys. Rep., 262, 1 

Sahni, V., & Starobinsky, A. 2006, Int. J. Mod. Phys. D, 15, 2105 

Samushia, L., & Ratra, B. 2006, ApJ, 650, L5 

Sanchez, A. G., et al. 2006, MNRAS, 366, 189 



-65- 



Sandage, A. 1978, AJ, 83, 904 

Sandage, A. R. 1995 in Sandage, A. R., Kron, R. G., k Longair, M. S., The Deep Universe: Saas-Fee 
Advanced Course 23, eds. B. Binggeli and R. Buser (Berlin: Springer), 1 

Sandage, A., et al. 2006, ApJ, 653, 843 

Sanders, R. H., k McGaugh, S. S. 2002, ARA&A, 40, 263 

Santiago, B. X., et al. 1995, ApJ, 446, 457 

Sasaki, S. 1996, PASJ, 48, L119 

Sato, K. 1981, Phys. Lett. B, 99, 66 

Sato, K. 1981, MNRAS, 195, 467 

Saunders, W., et al. 2000, MNRAS, 317, 55 

Schimd, C, et al. 2007, A&A, 463, 405 

Schneider, D. P., et al. 2005, AJ, 130, 367 

Schneider, P. 2006, in Kochanek, C. S., Schneider, P., k Wambsganss, J. 2006, 269 
Seldner, M., Siebers, B., Groth, E. J., k Peebles, P. J. E. 1977, AJ, 82, 249, 313, 314 
Seljak, U., Slosar, A., k McDonald, P. 2006, J. Cosmology Astropart. Phys., 0610, 014 



Sen, A. A., k Scherrer, R. J. 2007, astro-ph/0703416 



Shandarin, S. F., k Zel'dovich, Ya. B. 1989, Rev. Mod. Phys., 61, 185 
Shane, C. D., k Wirtanen, C. A. 1954, AJ, 59, 285 

Shapley, H., k Ames, A. 1932, Ann. Harvard Coll. Obs., 88, No. 2 (Cambridge: Astron. Obs. 
Harvard Coll.) 

Shectman, S. A., et al. 1996, ApJ, 470, 172 

Sheth, R. K., Mo, H. J., k Tormen, G. 2001, MNRAS, 323, 1 

Sheth, R. K., k Tormen, G. 1999, MNRAS, 308, 119 

Sheth, R. K., k van de Weygaert, R. 2004, MNRAS, 350, 517 

Silk, J. 1968, ApJ, 151, 459 

Silk, J. 1977, ApJ, 211, 638 



-66- 



Silk, J., & Rees, M. J. 1998, A&A, 331, LI 

Simon, J., Verde, L., & Jimenez, R. 2005, Phys. Rev. D, 71, 123001 
Slipher, V. M. 1917, Proc. Am. Philo. Soc, 56, 403 

Smirnov, Yu. N. 1964, English translation in 1965, Sov. Astron.-AJ, 8, 864 

Smith, R. W. 1982, The Expanding Universe: Astronomy's 'Great Debate' 1900-1931 (Cambridge: 
Cambridge Univ.) 

Smith, S. 1936, ApJ, 83, 23 

Smoot, G. F., et al. 1992, ApJ, 396, LI 

Somerville, R. S., & Primack, J. R. 1999, MNRAS, 310, 1087 
Spergel, D. N., et al. 2007, ApJS, 170, 377 

Spinrad, H. 2005, Galaxy Formation and Evolution (Berlin: Springer) 

Springel, V., Yoshida, N., & White, S. D. M. 2001, New A, 6, 79 

Springel, V., et al. 2005, New A, 435, 629 

Starobinsky, A. A. 1982, Phys. Lett. B, 117, 175 

Steidel, C. C, et al. 2004, ApJ, 604, 534 

Steigman, G. 2006, Int. J. Mod. Phys. E, 15, 1 

Steinhardt, P. J. 2003, Phil. Trans. R. Soc. Lond. A, 361, 2497 

Stoughton, C, et al. 2002, AJ, 123, 485 

Strauss, M. A., et al. 1992, ApJS, 83, 29 

Strauss, M. A., & Willick, J. A. 1995, Phys. Rep., 261, 271 

Subramanian, K. 2005, Curr. Sci., 88, 1068 

Subramanian, K. 2006, Astron. Nachr., 327, 403 

Sugiyama, N., & Gouda, N. 1992, Prog. Theor. Phys., 88, 803 

Swaters, R. A., Madore, B. F., & Trewhella, M. 2000, ApJ, 531, L107 

Sylos Labini, F., Montouri, M., & Pietronero, L. 1998 Phys. Rep., 293, 61 

Tasitsiomi, A. 2003, Int. J. Mod. Phys. D, 12, 1157 



-67- 



Tegmark, M., et al. 2004b, ApJ, 606, 702 

Tegmark, M., de Oliveira-Costa, A., & Hamilton, A. J. S. 2003, Phys. Rev. D, 68, 123523 

Tegmark, M., et al. 2006, Phys. Rev. D, 74, 123507 

Tegmark, M., et al. 1998, ApJ, 499, 555 

Tegmark, M., et al. 2004a, Phys. Rev. D, 69, 103501 

Tinker, J. L., Weinberg, D. H., & Zheng, Z. 2006, MNRAS, 368, 85 

Totsuji, H., & Kihara, T. 1969, PASJ, 21, 221 

Trimble, V. 1987, ARA&A, 25, 425 

Trodden, M. 2004, eConf C040802, L018 (Available at www.stanford.slac.edu/econf- 
/c040802/proceedings .htm.) 

Turner, E. L. 1990, ApJ, 365, L43 

Turner, M. S., & Widrow, L. M. 1988, Phys. Rev. D, 37, 2743 

Turok, N., ed. 1997, Critical Dialogues in Cosmology (Singapore: World Scientific) 

van de Hulst, H. C, Raimond, E., & van Woerden, H. 1957, Bull. Astron. Inst. Netherlands, 14, 1 

Verde, L., et al. 2002, MNRAS, 335, 432 

Verde, L., et al. 2003, ApJS, 148, 195 

Vettolani, G., et al. 1998, A&AS, 130, 323 

Villumsen, J. V. 1989, ApJS, 71, 407 

Vogeley, M. S., Park, C, Geller, M. J., & Huchra, J. P. 1992, ApJ, 391, L5 

Vogeley, M. S., et al. 1994, ApJ, 420, 525 

Vogeley, M. S., & Szalay, A. S. 1996, ApJ, 465, 34 

Voit, G. M. 2005, Rev. Mod. Phys., 77, 207 

Waga, I., & Frieman, J. A. 2000, Phys. Rev. D, 62, 043521 

Wagoner, R. V., Fowler, W. A., & Hoyle, F. 1967, ApJ, 148, 3 

Wake, D. A., et al. 2004, ApJ, 610, L85 

Wang, Y. 2006, ApJ, 647, 1 



-68- 



Wang, Y., & Mukherjee, P. 2006, ApJ, 650, 1 



Wang, Y., & Mukherjee, P. 2007, astro-ph/0703780 



Wang, Y., & Tegmark, M. 2005, Phys. Rev. D, 71, 103513 

Webb, S. 1999, Measuring the Universe: The Cosmological Distance Ladder (Chichester: Springer- 
Praxis) 

Weinberg, S. 1972, Gravitation and Cosmology: Principles and Applications of the General Theory 
of Relativity (New York: Wiley) 

Weinberg, S. 1993, The First Three Minutes: A Modern View of the Origin of the Universe, updated 
ed. (New York, Basic Books) 

Weyl, H. 1923, Phys. Zeit., 24, 230 (An English translation should appear in Gen. Relativ. Gravit. 
in 2007.) 

White, M., & Cohn, J. D. 2002, Am. J. Phys., 70, 106 
White, S. D. M. 1979, MNRAS, 186, 145 

White, S. D. M. 1996, in Cosmology and Large Scale Structure: Les Houches Session LX, eds. R. 
Schaeffer, J. Silk, M. Spiro, M., and J. Zinn-Justin (Amsterdam: Elsevier), 349 

White, S. D. M., Navarro, J. F., Evrard, A. E., & Frenk, C. S. 1993, Nature, 366, 429 

White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 

Widrow, L. M. 2002, Rev. Mod. Phys., 74, 775 

Wilkinson, D. T., & Peebles, P. J. E. 1990, in The Cosmic Microwave Background: 25 Years Later, 
eds. N. Mandolesi and N. Vittorio (Dordrecht: Kluwer), 17 

Willick, J. A., Strauss, M. A., Dekel, A., & Kolatt, T. 1997, ApJ, 486, 629 

Wilson, R. W. 1983, in Serendipitous Discoveries in Radio Astronomy, eds. K. Kellermann and B. 
Sheets (Green Bank: NRAO), 185 

Wilson, K. M., Chen, C, & Ratra, B. 2006, Mod. Phys. Lett. A, 21, 2197 

Wright, E. L. 2006, in Relativistic Astrophysics and Cosmology - Einstein's Legacy, ESO As- 
trophysics Symposia, eds. B. Aschenbach, V. Burwitz, G. Hasinger, and B. Leibundgut 
(Springer- Verlag), in press, [astro-ph/ 0603750 

Wright, E. L. 2007, ApJ, 664, 633 

Wright, E. L., et al. 1994, ApJ, 420, 450 



-69- 



Wu, J. H. P., et al. 2001, Phys. Rev. Lett., 87, 251303 
Younger, J. D., Bahcall, N. A., & Bode, P. 2005, ApJ, 622, 1 
Zee, A. 1980, Phys. Rev. Lett., 44, 703 
Zehavi, I., et al. 2005, ApJ, 630, 1 

Zel'dovich, Ya. B. 1963a, English translation in 1964, Sov. Phys. Usp., 6, 475 
Zel'dovich, Ya. B. 1963b, English translation in 1963, Atomic Energy, 14, 83 
Zel'dovich, Ya. B. 1965, Adv. A&A, 3, 241 
Zel'dovich, Ya. B. 1970, A&A, 5, 84 
Zel'dovich, Ya. B. 1972, MNRAS, 160, IP 

Zel'dovich, Ya. B., Kurt, V. G., & Sunyaev, R. A. 1968, English translation in 1969, Sov. Phys. 
JETP, 28, 146 

Zel'dovich, Ya. B., & Novikov, I. D. 1983, Relativistic Astrophysics. Vol. 2: The Structure and 
Evolution of the Universe (Chicago: Univ. of Chicago) 

Zhang, J., Zhang, X., & Liu, H. 2007, |astro-ph/0612642| 

Zhao, C, et al. 2007, Phys. Lett. B, 648, 8 

Zurek, W. H., Quinn, P. J., Salmon, J. K., & Warren, M. S. 1994, ApJ, 431, 559 

Zwicky, F. 1933, Helvetica Physica Acta, 6, 110 (An English translation should appear in Gen. 
Relativ. Gravit. in 2008.) 

Zwicky, F. 1937, ApJ, 86, 217 

Zwicky, F., Herzog, E., & Wild, P. 1961-1968, Catalogue of Galaxies and of Clusters of Galaxies 
(Pasadena: California Institute of Technology) 



This preprint was prepared with the AAS IATgjX macros v5.2. 



